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Abstract We review our algebraic framework for linear boundary problems (con- 
centrating on ordinary differential equations). Its starting point is an appropriate al- 
gebraization of the domain of functions, which we have named integro-differential 
algebras. The algebraic treatment of boundary problems brings up two new alge- 
braic structures whose symbolic representation and computational realization is 
based on canonical forms in certain commutative and noncommutative polyno- 
mial domains. The first of these, the ring of integro-differential operators, is used 
for both stating and solving linear boundary problems. The other structure, called 
integro-differential polynomials, is the key tool for describing extensions of integro- 
differential algebras. We use the canonical simplifier for integro-differential poly- 
nomials for generating an automated proof establishing a canonical simplifier for 
integro-differential operators. Our approach is fully implemented in the TH30REMV 
system; some code fragments and sample computations are included. 
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1 Introduction 

Overall View. When problems from Analysis — notably differential equations — are 
treated by methods from Symbolic Computation, one speaks of Symbolic Analysis, 
as in the eponymous workshops of the FoCM conference series [34]. Symbolic 
Analysis is based on algebraic structures, as all other symbolic branches, but its 
special flavor comes from its connection with analytic and numeric techniques. As 
most differential equations arising in the applications can only be solved numeri- 
cally, this connection is absolutely vital. 

If symbolic techniques cannot solve "most" differential equations, what else can 
they do? The answers are very diverse (reductions, normal forms, symmetry groups, 
singularity analysis, triangularization etc), and in the frame of this paper we can only 
point to surveys like [79] and [36, §2.11]. In fact, even the notion of "solving" is 
quite subtle and can be made precise in various ways. Often a symbolic method will 
not provide the "solution" in itself but valuable information about it to be exploited 
for subsequent numerical simulation. 

Our own approach takes a somewhat intermediate position while diverging rad- 
ically in another respect: Unlike most other symbolic methods known to us, we 
consider differential equations along with their boundary conditions. This is not 
only crucial for many applications, it is also advantageous from an algebraic point 
of view: It allows to define a linear operator, called the Green's operator, that maps 
the so-called forcing function on the right-hand side of an equation to the unique 
solution determined by the boundary conditions. This gives rise to an interesting 
structure on Green's operators and on boundary problems (Sect. 5). Algebraically, 
the consequence is that we have to generalize the common structure of differential 
algebras to what we have called integro-differential algebras (Sect. 3). 

Regarding the solvability issues, the advantage of this approach is that it uncou- 
ples the task of finding an algebraic representation of the Green's operator from that 
of carrying out the quadratures involved in applying the Green's operator to a forc- 
ing function. While the latter may be infeasible in a symbolic manner, the former 
can be done by our approach (with numerical quadratures for integrating forcing 
functions). 

The research program just outlined has been pursued in the course of the SFB 
project F013 (see below for a brief chronology), and the results have been reported 
elsewhere [70, 73, 66]. For the time being, we have restricted ourselves to linear 
boundary problems, but the structure of integro-differential polynomials [72] may 
be a first stepping stone towards nonlinear Green's operators. Since the algebraic 
machinery for Green's operators is very young, our strategy was to concentrate first 
on boundary problems for ordinary differential equations (ODEs), with some first 
steps towards partial differential equations (PDEs) undertaken more recently [74]. 
For an application of our methods in the context of actuarial mathematics, we refer 
to [2], for a more algebraic picture from the skew -polynomial perspective see [67]. 

New Results. In the present paper, we will present a new confluence proof for the 
central data structure used in our approach: As the algebraic language for Green's 
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operators, the integro-differential operators (Sect. 4) are defined as a ring of non- 
commutative polynomials in infinitely many variables, modulo an infinitely gener- 
ated ideal. While the indeterminates represent the basic operations of analysis (dif- 
ferentiation, integration, extraction of boundary values and multiplication by one of 
infinitely many coefficient functions), this ideal specifies their interaction (e.g. the 
fundamental theorem of calculus describing how differentiation and integration in- 
teract). Our new proof is fully automated within the TH30REMV system (Sect. 2), 
using a generic noncommutative polynomial reduction based on a noncommutative 
adaption of reduction rings [23]; see also [83] for a short outline of the proof. 

In a way, the new proof completes the circle started with the ad-hoc confluence 
proof in [69]. For the latter, no algebraic structure was available for coping with 
certain expressions that arise in the proof because they involved generic coefficient 
functions along with their integrals and derivatives (rather than the operator indeter- 
minates modeling integration and differentiation!), while this structure is now pro- 
vided by the afore-mentioned integro-differential polynomials (Sect. 6). Roughly 
speaking, this means within the spectrum between rewrite systems (completion by 
the Knuth-Bendix procedure) and Grobner bases (completion by Buchberger's algo- 
rithm), we have moved away from the former towards the latter [20]. We will come 
back to this point later (Sect. 7). 

Moreover, the paper includes the following improvements and innovations: The 
setting for Grobner bases and the Buchberger algorithm are introduced generically 
for commutative and noncommutative rings (allowing infinitely many variables and 
generators), based on reduction rings and implemented in the TH30REMV system 
(Sect. 2). The presentation of integro-differential algebras is streamlined and gener- 
alized (Sect. 3). For both of the main computational domains — integro-differential 
operators and integro-differential polynomials — we have a basis free description 
while a choice of basis is only need for deciding equality (Sects. 4, 6). The con- 
struction of integro-differential polynomials, which was sketched in [72], is carried 
out in detail (Sect. 6). In particular, a complete proof of the crucial result on canon- 
ical forms (Thm. 42) is now given. 

Chronological Outline. As indicated above, this paper may be seen as a kind of 
target line for the research that we have carried out within Project F1322 of the SFB 
F013 supported by the Austrian Science Fund (FWF). We have already pointed 
out the crucial role of analysis/numerics in providing the right inspirations for the 
workings of Symbolic Analysis. The development of this project is an illuminat- 
ing and pleasant case in point. It was initiated by the stimulating series of Hilbert 
Seminars conducted jointly by Bruno Buchberger and Heinz W. Engl from Octo- 
ber 2001 to July 2002, leading to the genesis of Project F1322 as a spin-off from 
Projects F1302 (Buchberger) and F1308 (Engl). Triggered by the paper [42], the 
idea of symbolic operator algebras emerged as a common leading theme. It en- 
gendered a vision of transplanting certain ideas like the Moore-Penrose inverse on 
Hilbert spaces from their homeground in functional analysis into a new domain 
within Symbolic Analysis, where powerful algebraic tools like Grobner bases are 
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available [17, 18, 9, 21]. This vision eventually crystallized in the algebraic machin- 
ery for computing Green's operators as described before. 

In the early stage of the project, those two main tools from analysis (Moore- 
Penrose inverse) and algebra (Grobner bases) were welded together in a rather ad- 
hoc manner, but it did provide a new tool for solving boundary problems [71]. In 
the course of the dissertation [69], a finer analysis led to a substantial simplification 
where the Moore-Penrose inverse was superseded by a purely algebraic formulation 
in terms of one-sided inverses and the expensive computation of a new noncom- 
mutative Grobner basis for each boundary problem was replaced by plain reduction 
modulo a fixed Grobner basis for modeling the essential operator relations. The re- 
sulting quotient algebra (called "Green's polynomials" at that time) is the precursor 
of the integro-differential operators described below (Sect. 4). The final step towards 
the current setup was the reformulation and generalization in a differential algebra 
setting [73] and in an abstract linear algebra setting [66]. 

The advances on the theoretical side were paralleled by an early implementa- 
tion of the algorithm for computing Green's operators. While the ad-hoc approach 
with computing Grobner bases per-problem was carried out by the help of NCAI- 
gebra, a dedicated Mathematica package for noncommutative algebra [42], the fixed 
Grobner basis for simplifying Green's operator was implemented in the TH30REMV 
system [26]; see Sect. 2 for a general outline of this system. As the new differential 
algebra setting emerged, however, it became necessary to supplant this implemen- 
tation by a new one. It was again integrated in the TH30REMV system, but now in 
a much more intimate sense: Instead of using a custom-tailored interface as in [69], 
the new package was coded directly in the TH30REMV language using the elegant 
structuring constructs of functors [25]. Since this language is also the object lan- 
guage of the provers, this accomplishes the old ideal of integrating computation and 
deduction. 

The presentation of several parts of this paper — notably Sects. 3, 4, 5 — benefited 
greatly from a lecture given in the academic year 2009/10 on Symbolic Integral 
Operators and Boundary Problems by the first two authors. The lecture was associ- 
ated with the Doctoral Program "Computational Mathematics: Numerical Analysis 
and Symbolic Computation" (W1214), which is a follow-up program to the SFB 
F013. We would like to thank our students for the lively discussions and valuable 
comments. 

Overview of the Paper. We commence by having a closer look at the TH30REMV 
system (Sect. 2), which will also be used in all sample computations presented in 
subsequent sections; both the sample computations and the TH30REMV program 
code is available in an executable Mathematica notebook from www . theorema . 
org. We discuss canonical simplifiers for quotient structures and Grobner bases in 
reduction rings, and we give a short overview of the functors used in building up 
the hierarchy of the algebraic structures used in the computations. The main struc- 
ture among these is that of an integro-differential algebra (Sect. 3), which is the 
starting point for the integro-differential operators as well as the integro-differential 
polynomials. Since the former are, in turn, the foundation for computing Green's op- 
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erators for boundary problems, we will next summarize the construction of integro- 
differential operators and their basic properties (Sect. 4), while the algorithms for 
solving and factoring boundary problems are explained and exemplified thereafter 
(Sect. 5). Driving towards the focus point of this paper, we describe then the alge- 
bra of integro-differential polynomials (Sect. 6), which will be the key tool to be 
employed for the confluence proof. Since this proof is reduced to a computation in 
TH30REMV, we will only explain the main philosophy and show some representa- 
tive fragments (Sect. 7). We wind up with some thoughts about open problems and 
future work (Sect. 8). 



2 Data Structures for Polynomials in Theorema 

The Theorema Functor Language. The TH30REMV system [26] was designed by 
B. Buchberger as an integrated environment for proving, solving and computing 
in various domains of mathematics. Implemented on top of Mathematica, its core 
language is a version of higher-order predicate logic that contains a natural pro- 
gramming language such that algorithms can be coded and verified in a unified 
formal frame. In this logic-internal programming language, functors are a power- 
ful tool for building up hierarchical domains in a modular and generic way. They 
were introduced and first implemented in TH30REMV by B. Buchberger. The gen- 
eral idea — and its use for structuring those domains in which Grobner bases can be 
computed — is described in [23, 25], where one can also find references to pertinent 
early papers by B. Buchberger. See also [87] for some implementation aspects of 
functor programming. 

The notion of functor in TH30REMV is akin to functors in ML, not to be con- 
fused with the functors of category theory. From a computational point of view, 
a TH30REMV functor is a higher-order function that produces a new domain from 
given domains, where each domain is considered as a bundle of operations (includ- 
ing relations qua boolean-valued operations — in particular also carrier predicates). 
Operations in the new domain are defined in terms of operations in the underlying 
domains. 

Apart from this computational aspect, functors also have an important reason- 
ing aspect — a functor transports properties of the input domains to properties of the 
output domain, typical examples being the various "preservation theorems" in math- 
ematics: "If R is a ring, then R[x] is also a ring". This means the functor R i-)- R[x] 
preserves the property of being a ring, in other words: it goes from the "category of 
rings" to itself. In this context, a category is simply a collection of domains charac- 
terized by a common property (a higher-order predicate on domains). 

See below for an example of a functor named LexWords. It takes a linearly or- 
dered alphabet L as input domain and builds the word monoid over this alphabet: 
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Definition^ "Word Monoid", any[L], 
LexWords [L] =Functor^W, any[v, w, £, rj, £, 77], 

8 = <> 



If is-tuple[w] 
A L, v ,M? [wj 



n = <> 

w 

V * W = V B W 

w 

[<ri, j)> > <>) ~ True 
(<> > <f7>) o False 



]] 



V / 



f) A <*»> > <^> 



Here %, fj are sequence variables, i.e. they can be instantiated with finite sequences 
of terms. The new domain W has the following operations: W[g] denotes the carrier 
predicate, the neutral element is given by W[d], the multiplication W[*] is defined as 
concatenation, and W[>] defines the lexicographic ordering on W. 

In the following code fragments, we illustrate one way of building up polyno- 
mials in TH30REMV starting from the base categories of fields with ordering and 
ordered monoids. Via the functor FreeModule, we construct first the free vector 
space V over a field K generated by the set of words in an ordered monoid W. The 
elements of V are described by V[g] as lists of pairs, each pair containing one (non- 
zero) coefficient from K and one basis vector from W, where the basis vectors are 
ordered according to the ordering on W. The operations of addition, subtraction and 
scalar multiplication are defined recursively, using the operations on K and W: 

Definition ^ "Free Module", any[K, W] , 
FreeModule [K, W] = Functor^V, any[c, d, x, y, r], A, it, y] , 

s = <> 



6[X] < 



where = \x\ , f\ 



is-tuple [x] 



A 



is-tupletx^ 

1 *i 1 = 2 

is-coef f [ (x ± ) 



is-bvec[ (x L ) 2 ] 

v 

(Xi) 2 > (Xi.,i) 2 



is-bvec[ff] o € [{] 
v w 

is-coef f [c] o |e [c] A c * °) 
= <> 



<<c, £>, x> + <<d, r,), y> 



<c, £>- (<x> j <<d, n>, y>) 
<d, r,) - (<<c, £>, x> + <y>) 

(ojd. e)-(<x> ; <y>) 
<x> + <y> 



H 

n>£ 

W 

(f = '7)A c + d ' !0 



<= otherwise 
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x • = <> 

V 

]] 

<<d, /;>, y> = ( c * d / i) <y> 
<<c, O, *>^d= (c*d, f)-<x> -d 

By the MonoidAlgebra functor we extend this domain, introducing a multipli- 
cation using the corresponding operations in K and W: 

MonoidAlgebra[K, W] = where = FreeModule [K, W] , 
Functor^P, any[c, d, f, g, £, J7, in, n] , 

s = <) 

...(* linear operations from V *) 
(* multiplication *) 

<>*g= <> 

P 

f *<> = <> 

P 

<<c, £>, fi>.<<d, ,7), S> = ((c.d, £.17)) + <<c, f>> »<S> + <i> *<<d, r,>, S> 

P \\K WMP PPP 

The new domain inherits the structure on the elements of V. 

The main advantage of the above construction is that it is fully generic: Not only 
can it be instantiated for different coefficient rings (or fields) and different sets of 
indeterminates, it comprises also the commutative and noncommutative case (where 
W is instantiated respectively by a commutative and noncommutative monoid). 

Quotient Structures and Canonical Simplifiers. In algebra (and also in the rest of 
mathematics), one encounters quotient structures on many occasions. The general 
setting is a set A with various operations (an algebra in the general sense used in 
Sect. 6) and a congruence relation = on A, meaning an equivalence relation that 
is compatible with all the operations on A. Then one may form the quotient A/=, 
which will typically inherit some properties of A. For example, A/= belongs to the 
category of rings if A does, so we can view the quotient construction A h-> A/= as a 
functor on the category of rings. 

But for computational purposes, the usual set-theoretic description of A/= as a 
set of equivalence classes is not suitable (since each such class is typically uncount- 
ably infinite). We will therefore use an alternative approach that was introduced 
in [27] as a general framework for symbolic representations. The starting point is a 
canonical simplifier forA/=, meaning a map a : A — s- A such that 

a(a)=a and o(a) = o(a') whenever a = a 1 . (1) 

The set A = cr(A) is called the associated system of canonical forms for A/=. 

Clearly canonical simplifiers exist for every quotient A/=, but for computational 
purposes the crucial question is whether a is algorithmic. Depending on A/=, it 
may be easy or difficult or even impossible to construct a computable a : A — 5- A. In 
the examples that we will treat, canonical simplifiers are indeed available. 



]] 
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Canonical simplifiers are also important because they allow us to compute in the 
quotient structure. More precisely, one can transplant the operations on A to A by 
defining (0{a\, . . . ,a n ) = a(co(a\,. . . ,a n )) for every operation co on A. With these 
new operations, one may easily see that A is isomorphic to the quotient A/=; see 
the Theorem "Canonical simplification and computation" in [27, p. 13]. 

There is an intimate relation between canonical forms and normal forms for 
rewrite systems (Sect. 4 contains some basic terminology and references). In fact, 
every rewrite system — > on an algebraic structure A creates an equivalence relation 
=, the symmetric closure of A. Thus a = a' if and only if a and a' can be connected 
by an equational chain (using the rewrite rules in either direction). Typically, the 
relation = will actually be a congruence on A, so that the quotient A/= has a well- 
defined algebraic structure. Provided the rewrite system is noetherian, the normal 
forms of — > are then also canonical forms for A/=. Hence we will often identify 
these terms in a rewriting context. 

For our implementation, we use canonical simplifiers extensively. In fact, the 
observation made above about computing in the quotient structure is realized by a 
TH30REMV functor, which is applied at various different places. Here A is typically 
a /f-algebra, with the ground field K being Q or computable subfields of R and (D. 

Reduction Rings and Grobner Bases. For defining reduction on polynomials, we 
use the reduction ring approach in the sense of [19, 23]. For commutative reduction 
rings, see also [81, 82]; for another noncommutative approach we refer to [57, 58, 



To put it simply, a reduction ring is a ring in which Grobner bases can be done. 
A full axiomatization for the commutative case is given in [19]. If such rings satisfy 
certain additional axioms (defining the category of so-called "Grobner rings"), then 
Grobner bases can be computed by iterated S-polynomial reduction in the given 
ring — this is the Grobner Ring Extension Theorem, stated and proved in [19]. 

A detailed presentation of their construction in the TH30REMV setting was given 
in [22, 24]; it is the starting point for our current work. At this point we do not give 
an axiomatic characterization for noncommutative reduction rings, but we do use 
a construction that is similar to the commutative setting. Thus we endow a poly- 
nomial domain P, built via the MonoidAlgebra functor with word monoid W and 
field K, with the following three operations: a noetherian (partial) ordering, a binary 
operation least common reducible, and a binary operation reduction multiplier. The 
noetherian ordering is defined in the usual way in terms of the given orderings on K 
and W. 

The basic idea of reduction multipliers is to answer the question: "With which 
monomial do I have to multiply a given polynomial so that it cancels the leading 
term of another given polynomial?" In the noncommutative case, the corresponding 
operation rdm splits into left reduction multiplier lrdm and its right counterpart 
rrdm defined as follows: 



59]. 
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rrdm[<<c, f > , m> , <<d, 77), n>] = 




K w / / K K 



P 







otherwise 



p 



Here the divisibility relation | on W checks whether a given word occurs within 
another word, and the corresponding quotients lquot and rquot yield the word 
segments respectively to the left and to the right of this occurrence. Since the scalars 
from K commute with the words, it is an arbitrary decision whether one includes it 
in the right (as here) or left reduction multiplier. In typical cases, this scalar factor 
is just rdm[c,d] = c/d. 

The operations relating Grobner bases are introduced via a functor which is 
called GroebnerExtension. It defines polynomial reduction using reduction mul- 
tipliers (note that this includes also the commutative case, where one actually needs 
only one reduction multiplier, the other one being unity): 



The next step is to introduce reduction modulo a system of polynomials. For some 
applications (like the integro-differential operators described in Sect. 4), it is neces- 
sary to deal with infinite reduction systems: the polynomial ring contains infinitely 
many indeterminates, and reduction is applied modulo an infinite set of polynomi- 
als. In other words, we want to deal with an infinitely generated ideal in an infinitely 
generated algebra. 

This is a broad topic, and we cannot hope to cover it in the present scope. In gen- 
eral one must distinguish situations where both the generators of the ideal and the 
algebra are parametrized by finitely many families involving finitely many param- 
eters and more general algebras/ideals where this is not so. In the latter case, one 
works with finite subsets, and all computations are approximate: one never catches 
the whole algebraic picture. Fortunately, the applications we have in mind — in par- 
ticular the integro-differential operators — are of the first type where full algorithmic 
control can be achieved. However, most of the common packages implementing 
noncommutative Grobner bases do not support such cases [55, 56]. For some recent 
advances, we refer the reader to [3, 14, 43, 51] as well as Ufnarovski's extensive 
survey chapter [86]. 

Let us point out just one important class of decidable reductions in infinitely 
generated algebras — if an infinite set of (positively weighted) homogeneous poly- 
nomials is given, which is known to be complete for each given degree (see [51] 
for the proof) since one can compute a truncated Grobner basis of such a graded 
ideal, which is finite up to a given degree. But if the given set is not homogeneous 
or cannot be clearly presented degree by degree, basically nothing can be claimed 
in general. Unfortunately, the applications we have in mind seem to be of this type. 

In our setting, infinitely generated ideals are handled by an algorithmic operation 
for instantiating reduction rules. The reduction of polynomial f modulo a system S 
is realized thus: 



hred[f, g] =f-lrdm[f, g]*g*rrdm[f, g] 

G p P P P P 



G p P P P P 



hredp [f, 1, g, r] =f-l*g*r 

hred[f, S] = where q = S [f] , hredp[f, qi, q2, q 3 ] 
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where S[f ] is the operation that decides if there exists g modulo which f can be re- 
duced, and it returns a triple containing the g and the left/right reduction multipliers 
needed for performing the reduction. 

The main tool for the Grobner bases construction, namely the notion of S- 
polynomial, can now be defined in terms of the least common reducible: 

spol[f, g] = where L = lcrd[f, g] , hredp[L, f] -hredp[L, g] 

G L p G " a 1 

Here lcrd[f ,g] represents the smallest monomial that can be reduced both mod- 
ulo f and modulo g, built from the least common reducible of the corresponding 
coefficients in K and the least common multiple of the words in W: 

lcrd[<<c, £>, fi>, <<d, n >, n>] = ^lcrd[c, d] , lcm[f, 77] 

In our setting, the lcrd[c, d] can of course be chosen as unity since we work over a 
field K, but in rings like TL one would have to use the least common multiple. 

Finally, Grobner bases are computed by the usual accumulation of S-polynomials 
reduction, via the following version of Buchberger algorithm [17]: 

Gb[R, S] = where [pairs = /{Ri, Rj) | RijlRj), Gb[R, pairs, S] 1 

G > i=l,..,|Rl ' G 

j.l,.., |R| 

Gb[R, <>, S] = R 

G 

Gb [R, <<f, g>, m>, S] = where [h = tredjspol [f, g] , s] , 
Gb[R, <m>, S] •= h = 

G P 

Gb[R-h, (|h-»S) x <m>) = (R «- h) , S] <= otherwise 

G 
] 

Total reduction modulo a system, denoted here by tred, is computed by iteratively 
performing reductions, until no more reduction is possible. The above implementa- 
tion of Buchberger's algorithm is again generic since it can be used in both commu- 
tative and noncommutative settings. For finitely many indeterminates, the algorithm 
always terminates in the commutative case (by Dickson's Lemma); in the noncom- 
mutative setting, this cannot be guaranteed in general. For our applications we also 
have to be careful to ensure that the reduction systems we use are indeed noetherian 
(Sect. 4). 



3 Integro-Differential Algebras 

For working with boundary problems in a symbolic way, we first need an algebraic 
structure having differentiation along with integration. In the following definitions, 
one may think of our standard example & = C°° (R), where d = ' is the usual deriva- 
tion and J the integral operator 



/-> Caws 
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for a fixed a € R. 

Axioms and Basic Properties. Let A" be a commutative ring. We first recall that 
(^,d) is a differential K-algebra if d : & — > & is a /T-linear map satisfying the 
Leibniz rule 

d(fg)=d(f)g + fd(g). (2) 

For convenience, we may assume K < and we write /' as a shorthand for <?(/). 
The following definition [73] captures the algebraic properties of the Fundamental 
Theorem of Calculus and Integration by Parts. 

Definition 1. We call (J?,<3, J) an integro-differential algebra if (i^,<3) is a com- 
mutative differential A-algebra and J is a ^-linear section (right inverse) of d, i.e. 

(//)' = /, 0) 
such that the differential Baxter axiom 

Uf')(h')+I(fgy = (If')g+f(k') (4) 

holds. 

We refer to d and J respectively as the derivation and integral of and to (3) as 
section axiom. Moreover, we call a section J of d an integral for d if it satisfies (4). 
For the similar notion of differential Rota-Baxter algebras, we refer to [39] but see 
also below. 

Note that we have applied operator notation for the integral; otherwise, for ex- 
ample, the section axiom (3) would read (/(/))' = /, which is quite unusual at least 
for an analyst. We will likewise often use operator notation for the derivation, so 
the Leibniz rule (2) can also be written as dfg = (df)g + f(dg). For the future we 
also introduce the following convention for saving parentheses: Multiplication has 
precedence over integration, so jfjg is to be parsed as J(fjg). 

Let us also remark that Definition 1 can be generalized: First, no changes are 
needed for the noncommutative case (meaning J£" is noncommutative). This would 
for example be an appropriate setting for matrices with entries in = C°°[a,b], 
providing an algebraic framework for the results on linear systems of ODEs. Second, 
one may add a nonzero weight in the Leibniz axiom, thus incorporating also discrete 
models where d is the difference operator defined by {df)k = fk+i —fk- The nice 
thing is that all other axioms remain unchanged. For both generalizations confer 
also to [39]. 

We study first some direct consequences of the section axiom (3). For further 
details on linear left and right inverses, we refer for example to [ 1 3, p. 2 1 1 ] or to [63] 
in the context of generalized inverses. We also introduce the following names for the 
projectors and modules associated with a section of a derivation. 

Definition 2. Let (J^,<9) be a differential A-algebra and J a ^-linear section of d. 
Then we call the projectors 
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J = J o d and e = 1 — J o d 

respectively the initialization and the evaluation of & ', Moreover, we refer to 

^ = Ker(<9) = Ker(j) =Im(E) and J = Im(J) = Im( j) = Ker(E) 

as the submodules of respectively constant and initialized functions. 

Note that they are indeed projectors since joj = J'o(<9oJ')o<9 = jby (3), 
which implies eoe = 1- j- j+ joj = e. Asis well known [13, p. 209], every 
projector is characterized by its kernel and image — they form a direct decomposition 
of the module into two submodules, and every such decomposition corresponds to 
a unique projector. We have therefore a canonical decomposition 

& = tf + J r , 

which allows to split off the "constant part" of every "function" in & ' . 

Before turning to the other axioms, let us check what all this means in the stan- 
dard example & = C°°(R) with d = ^ and J = J*. Obviously, the elements of < io 
are then indeed the constant functions f(x) = c, while J? consists of those functions 
that satisfy the homogeneous initial condition f(a) = 0. This also explains the ter- 
minology for the projectors: Here e/ = f(a) evaluates / at the initialization point 
a, and jf = f — f(a) enforces the initial condition. Note that in this example the 
evaluation e is multiplicative; we will show below that this holds in any integro- 
differential algebra. 

The Leibniz rule (2) and the differential Baxter axiom (4) entail interesting prop- 
erties of the two submodules ^ and J^. For understanding these, it is more economic 
to forget for a moment about integro-differential algebras and turn to the following 
general observation about projectors on an algebra. We use again operator notation, 
giving precedence to multiplication over the linear operators. 

Lemma 3. Let E and J be projectors on a K-algebra with E +J = 1, set 

C = ha(E) = Ker(7) and I = Ker(£) = Jm(J). 

Then the following statements are equivalent: 

1. The projector E is multiplicative, meaning Efg = (Ef)(Eg). 

2. The projector J satisfies the identity (J f) (Jg) +Jfg= (Jf)g + f(Jg) ■ 

3. The submodule C is a subalgebra and the submodule I an ideal. 

Proof. L_Q_2. Multiplicativity of E = 1 — / just means 

fg - Jfg = fg - (Jf)g - f(Jg) + W) (Jg) • 

1. => 3. This follows immediately because C is the image and / the kernel of the 
algebra endomorphism E. 
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3. => 1. Let /, g be arbitrary. Since the given /^-algebra is a direct sum of C and /, we 
have / = f c +fi and g = g c + gi for f c = Ef,g c = EgeC and // = Jf, gI =Jg£l. 
Then 

Efg = Efcgc + Efcgi + Efigc + Efigi 

Since / is an ideal, the last three summands vanish. Furthermore, C is a subalgebra, 
so fcgc G C. This implies Efcgc = fcgc because £ is a projector onto C. □ 

This lemma is obviously applicable to integro-differential algebras & with the 
projectors E = e and J = j and with the submodules C = c & and I = J? because the 
differential Baxter axiom (4) is exactly condition 2. From now on, we will therefore 
refer to c € as the algebra of constant functions and to ,f as the ideal of initialized 
functions. Moreover, we note that in any integro-differential algebra the evaluation 
e = 1 — J o <9 is multiplicative, meaning 

E/g =(«/)(■*). (5) 

Altogether we obtain now the following characterization of integrals (note that the 
requirement that c € be a subalgebra already follows from the Leibniz axiom). 

Corollary 4. Let {^P ,d) be a differential algebra. Then a section J of d is an inte- 
gral if and only if its evaluation e = 1 — J o d is multiplicative, and if and only if 
,y = Im(J) is an ideal. 

Note that the ideal J? corresponding to an integral is in general not a differential 
ideal of We can see this already in the standard example C°°[0, 1], where J* 
consists of all / G C°°[0, 1] with /(0) = 0. Obviously J* is not differentially closed 
since xg / but x' = 1 £ J . 

The above corollary implies immediately that an integro-differential algebra & 
can never be a field since then the only possibilities for J would be and & . 
The former case is excluded since it means that Ker(<9) = J£\ contradicting the 
surjectivity of d. The latter case corresponds to Ker(<9) = 0, which is not possible 
because d 1 = 0. 

Corollary 5. An integro-differential algebra is never afield. 

In some sense, this observation ensures that all integro-differential algebras are 
fairly complicated. The next result points in the same direction, excluding finite- 
dimensional algebras. 

Proposition 6. The iterated integrals l,Jl,JjT,... are all linearly independent 
over K. In particular, every integro-differential algebra is infinite- dimensional. 

Proof. Let (u„) be the sequence of iterated integrals of 1. We prove by induction 
on n that UQ,Ui,...,u n are linearly independent. The base case n = is trivial. For 

the induction step from n ton+ 1, assume cquq^ \-c n +\u n +\ = 0. Applying <9" +1 

yields c n +\ = 0. But by the induction hypothesis, we have then also cq = ■ ■ ■ = c„ = 
0. Hence mq, ■ ■ ■ , « n +i ar e linearly independent. □ 
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Let us now return to our discussion of the differential Baxter axiom (4). We will 
offer an equivalent description that is closer to analysis. It is more compact but 
less symmetric. (In the noncommutative case one has to add the opposite version — 
reversing all products — for obtaining equivalence.) 

Proposition 7. The differential Baxter axiom (4) is equivalent to 



in the presence of the Leibniz axiom (2) and the section axiom (3). 

Proof. For proving (6) note that since J? is an ideal, fjg is invariant under the 
projector j and thus equal to j(fjg)' = jf'jg + jfg by the Leibniz axiom (2) 
and the section axiom (3). Alternatively, one can also obtain (6) from (4) if one 
replaces g by jg in (4). Conversely, assuming (6) we see that J? is an ideal of 
so Corollary 4 implies that J satisfies the differential Baxter axiom (4). □ 

For obvious reasons, we refer to (6) as integration by parts. The usual formula- 
tion J/G' = fG — If'G is only satisfied "up to a constant", or if one restricts G to 
Im(J). Substituting G = jg then leads to (6). But note that we have now a more 
algebraic perspective on this well-known identity of Calculus: It tells us how J" is 
realized as an ideal of & . 

Sometimes a variation of (6) is useful. Applying J to the Leibniz axiom (2) and 
using the fact that e = 1 — j is multiplicative (5), we obtain 



which we call the evaluation variant of integration by parts (a form that is also used 
in Calculus). Observe that, we regain integration by parts (6) upon replacing g by 
Jg in (7) since e jg = 0. 

Note that in general one cannot extend a given differential algebra to an integro- 
differential algebra since the latter requires a surjective derivation. For example, in 
(K[x 2 ],xd) the image of d does not contain 1. As another example (cf. Sect. 6), 
the algebra of differential polynomials & = K{u} does not admit an integral in the 
sense of Definition 1 since the image of d does not contain u. 

How can we isolate the integro part of an integro-differential algebra? The disad- 
vantage (and also advantage!) of the differential Baxter axiom (4) is that it entangles 
derivation and integral. So how can one express "integration by parts" without re- 
ferring to the derivation? 

Definition 8. Let & be a /f-algebra and J a /^-linear operation satisfying 



fjg=if8+jrjg, 



(6) 



fftf=fg-Sfg-{*f)(*g), 



(7) 



Uf)(h) = IfJg+JgJf- 



(8) 



Then (j£ , J) is called a Rota-Baxter algebra (of weight zero). 



Rota-Baxter algebras are named after Glen Baxter [7] and Gian-Carlo Rota [75]; 
see also [37, 38] for further details. In the following, we refer to (8) as Baxter axiom; 
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in contrast to the differential Baxter axiom (4), we will sometimes also call it the 
pure Baxter axiom. 

One might now think that an integro-differential algebra (J£",<9, J) is a differ- 
ential algebra (J^,<5) combined with a Rota-Baxter algebra (J^", J) such that the 
section axiom (3) is satisfied. In fact, such a structure was introduced independently 
by Guo and Keigher [39] under the name differential Rota-Baxter algebras. But we 
will see that an integro-differential algebra is a little bit more — this is why we also 
refer to (8) as "weak Baxter axiom" and to (4) and (6) as "strong Baxter axioms". 

Proposition 9. Let (&,d) be a differential algebra and J a section for d. Then J 
satisfies the pure Baxter axiom (8) if and only if J? = Im(J) is a subalgebra of ' & ' . 
In particular, j) is a Rota-Baxter algebra for any integro-differential algebra 
(&,dj). 

Proof. Clearly (8) implies that J* is a subalgebra of Conversely, if (jf)(jg) is 
contained in J , it is invariant under the projector j and must therefore be equal to 
Sd(Sf) (jg) = Sfjg + SgSf by the Leibniz axiom (2). □ 

So the strong Baxter axiom (4) requires that J* be an ideal, the weak Baxter 
axiom (8) only that it be a subalgebra. We will soon give a counterexample for 
making sure that (4) is indeed asking for more than (8), see Example 14. But before 
this we want to express the difference between the two axioms in terms of a linearity 
property. Recall that both d and J were introduced as TT-linear operations on & . 
Using the Leibniz axiom (2), one sees immediately that d is even ^-linear. It is 
natural to expect the same from J, but this is exactly the difference between (4) 
and (8). 

Proposition 10. Let ,d) be a differential algebra and J a section for d. Then j 
satisfies the differential Baxter axiom (4) if and only if it satisfies the pure Baxter 
axiom (8) and is c & -linear. 

Proof. Assume first that J satisfies the differential Baxter axiom (4). Then the pure 
Baxter axiom (8) holds by Proposition 9. For proving jcg = c jg for all and 
g £ J?, we use the integration-by-parts formula (6) and c' = 0. 

Conversely, assume the pure Baxter axiom (8) is satisfied and J is ^-linear. By 
Proposition 7 it suffices to prove the integration-by-parts formula (6) for f,g£ & . 
Since ^ = ^+/,we may first consider the case / £ c £ and then the case f & JP. 
But the first case follows from c £ -linearity; the second case means /= J/ for / G 
and (6) becomes the pure Baxter axiom (8) for / and g. □ 

Let us now look at some natural examples of integro-differential algebras, in 
addition to our standard examples C°°(R) and C°°[a,b]. 

Example 11. The analytic functions on the real interval [a,b] form an integro- 
differential subalgebra C m [a,b] of C°°[a,£>] over K = 1R or K = <D. It contains in 
turn the integro-differential algebra K[x,e Kx ] of exponential polynomials, defined as 
the space of all /f-linear combinations of x n e^ x , with n£W and A £ K. Finally, 
the algebra of ordinary polynomials K[x] is an integro-differential subalgebra in all 
cases. □ 
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All the three examples above have algebraic analogs, with integro-differential 
structures defined in the expected way. 

Example 12. For a field K of characteristic zero, the formal power series K[[x]] 
are an integro-differential algebra. One sets dx k = kx k ~ 1 and §x k = x k+l / (k + 1); 
note that the latter needs characteristic zero. The formal power series contain 
a highly interesting and important integro-differential subalgebra: the holonomic 
power series, defined as those whose derivatives span a finite-dimensional /f-vector 
space [29, 77]. 

Of course K[[x]] also contains (an isomorphic copy of) the integro-differential 
algebra of exponential polynomials. In fact, one can define K[x,e Kx ] algebraically 
as a quotient of the free algebra generated by the symbols x k and e^ x , with X rang- 
ing over K. Derivation and integration are then defined in the obvious way. The 
exponential polynomials contain the polynomial ring K[x] as an integro-differential 
subalgebra. When K — M, or K — (D, we use the notation K[x] and K[x,e ] both for 
the analytic and the algebraic object since they are isomorphic. □ 

The following example is a clever way of transferring the previous example to 
coefficient fields of positive characteristic. 

Example 13. Let K be an arbitrary field (having zero or positive characteristic). 
Then the algebra H(K) of Hurwitz series [46] over K is defined as the /f-vector 
space of infinite /f-sequences with the multiplication defined as 

(a„) ■ (b„) = (jl a i b n-ij 

for all (a„),(b n ) G H(K). If one introduces derivation and integration by 

d(a ,ai,a 2 ,...) = (a h a 2 ,...), 
J(flo,fli,---) = (0,a o ,ai,...), 

the Hurwitz series form an integro-differential algebra (H(K),d, J), as explained 
by [47] and [37]. Note that as an additive group, H(K) coincides with the formal 
power series ^[[z]], but its multiplicative structure differs: We have an isomorphism 

52 a„z n H> (n\a n ) 

from K[[z]] to H(K) if and only if K has characteristic zero. The point is that one 
can integrate every element of H(K), whereas the formal power series z p ~ l does not 
have an antiderivative in K[[z]] if K has characteristic p > 0. □ 

Now for the promised counterexample to the claim that the section axiom would 
suffice for merging a differential algebra {&,d) and a Rota-Baxter algebra (J£", J) 
into an integro-differential algebra (J^",<9, J). 
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Example 14. Set R = K\y]/(y 4 ) for K a field of characteristic zero and define d on 
& = R[x] as usual. Then (J^,<9) is a differential algebra. Let us define a /f-linear 
map J on ^ by 

J/ = f/+/(0,0)y 2 , (9) 

where J" is the usual integral on /?[jc] with.* 1 i— >x k+l /{k-\- 1). Since the second term 
vanishes under d, we see immediately that J is a section of d. For verifying the pure 
Baxter axiom (8), we compute 

(J/)(J» = (f/)(fs) +J 2 f («(0,0)/ + /(0,0)g) + /(0,0)g(0,0)y 4 , 
J/J^ - J/(jV + ^(0,0)y 2 ) = jVj«+*(0,0)y 2 jV- 

Since y 4 = and the ordinary integral J* fulfills the pure Baxter axiom (8), this 
implies immediately that J does also. However, it does not fulfill the differential 
Baxter axiom (4) because it is not c € -linear: Observe that c <§ is here Ker(<9) = R, so 
in particular we should have J(y - 1 ) = y ■ Jl . But one checks immediately that the 
left-hand side yields xy, while the right-hand side yields xy + y 3 . □ 

Ordinary Integro-Differential Algebras. The following example shows that our cur- 
rent notion of integro-differential algebra includes also algebras of "multivariate 
functions". 

Example 15. Consider & = C°° (R 2 ) with the derivation du = u x + u y . Finding sec- 
tions for d means solving the partial differential equation u x + u y = f. Its general 
solution is given by 

u( x ,y)= [ f(t,y-x+t)dt+g(y-x), 

J a 

where g G C°°(1R) and a € R are arbitrary. Let us choose a = for simplicity. In 
order to ensure a linear section, one has to choose g = 0, arriving at 

J/ = / f(t,y-x+t)dt, 

Using a change of variables, one may verify that J satisfies the pure Baxter ax- 
iom (8), so (J£", J) is a Rota-Baxter algebra. 

We see that the constant functions ^ = Ker(<3) are given by (x,y) i-tc(x—y) with 
arbitrary c G C°°(]R), while the initialized functions = Im(J) are those F G -!P that 
satisfy F(Q,y) = for all y G R. In other words, c <§ consists of all functions constant 
on the characteristic lines x— y = const, and of those satisfying the homogeneous 
initial condition on the vertical axis (which plays the role of a "noncharacteristic 
initial manifold"). This is to be expected since J integrates along the characteristic 
lines starting from the initial manifold. The evaluation e : J? — 5- & maps a function 
/ to the function (x,y) n> f(0,y— x). This means that / is "sampled" only on the 
initial manifold, effectively becoming a univariate function: the general point (x,y) 
is projected along the characteristics to the initial point (0,y — x). 
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Since e is multiplicative on Lemma 3 tells us that (jF,<9,J) is in fact an 
integro-differential algebra. Alternatively, note that J 1 is an ideal and that J is c €- 
linear. Furthermore, we observe that here the polynomials are given by K[x]. □ 

In the following, we want to restrict ourselves to boundary problems for ordinary 
differential equations. Hence we want to rule out cases like Example 15. The most 
natural way for distinguishing ordinary from partial differential operators is to look 
at their kernels since only the former have finite-dimensional ones. Note that in 
the following definition we deviate from the standard terminology in differential 
algebra [48, p. 58], where ordinary only refers to having a single derivation. 

From now on, we restrict the ground ring K to afield. We can now characterize 
when a differential algebra is ordinary by requiring that c £ be one-dimensional over 
K, meaning = K. 

Definition 16. A differential algebra (J^,<3) is called ordinary if dim^Ker(5) = 1. 

Note that except for Example 15 all our examples have been ordinary integro- 
differential algebras. The requirement of ordinariness has a number of pleasant 
consequences. First of all, the somewhat tedious distinction between the weak and 
strong Baxter axioms disappears since now & is an algebra over its own field of 
constants K = . Hence J is by definition "^-linear, and Lemma 10 ensures that 
the pure Baxer axiom (8) is equivalent to the differential Baxter axiom (4). Let us 
summarize this. 

Corollary 17. In an ordinary integro-differential algebra, the constant functions co- 
incide with the ground field, and the strong and weak Baxter axioms are equivalent. 

Recall that a character on an algebra (or group) is a multiplicative linear func- 
tional; this may be seen as a special case of the notion of character in representation 
theory, namely the case when the representation is one-dimensional. In our context, 
a character on an integro-differential algebra is a Zf-linear map <p : — > K sat- 
isfying (p(fg) = <p(f) q>(g) and a fortiori also <p(l) = 1. So we just require <p to be 
a /if-algebra homomorphism, as for example in [52, p. 407]. 

Ordinary integro-differential algebras will always have at least one character, 
namely the evaluation: One knows from Linear Algebra that a projector P onto 
a one-dimensional subspace [w] of a ZT-vector space V can be written as P(v) = 
(p(v)w, where <p : V —> K is the unique functional with (p(w) = 1 . If V is moreover a 
/if-algebra, a projector onto K = [1] is canonically described by the functional q> with 
normalization <p(l) = 1. Hence multiplicative projectors like e can be viewed as 
characters. In the next section, we consider other characters on for the moment 
let us note e as a distinguished character. We write & % for the set of all nonzero 
characters on a Zf-algebra JF, in other words all algebra homomorphisms & — > K. 

One calls a /f-algebra augmented if there exists a character on it. Its kernel J? 
is then known as an augmentation ideal and forms a direct summand of K\ see for 
example [32, p. 132]. Augmentation ideals are always maximal ideals (generalizing 
the C°°[a,b) case) since the direct sum J? = K + J 1 induces a ring isomorphism 
,^/,^ = K. Corollary 4 immediately translates to the following characterization of 
integrals in ordinary differential algebras. 
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Corollary 18. In an ordinary differential algebra (j?,<9), a section J of d is an 
integral if and only if its evaluation is a character if and only if J' = Im(J) is an 
augmentation ideal. 



Initial Value Problems. It is clear that in general we cannot assume that the solutions 
of a differential equation with coefficients in & are again in & . For example, in 
& = K[x], the differential equation u' — u = has no solution. In fact, its "actual" 
solution space is spanned by u(x) = e x if K — Ror K = (D. So in this case we should 
have taken the exponential polynomials & = K[x,e Kx ] for ensuring that u £ & . 
But if this is the case, we can also solve the inhomogeneous differential equation 
u' — u = f whose general solution is Ke x + e*Je~ x f, with J = j* as usual. Of course 
we can also incorporate the initial condition m(0) = 0, which leads to u = e?[e~ x f. 

This observation generalizes: Whenever we can solve the homogeneous differ- 
ential equation within J?, we can also solve the initial value problem for the corre- 
sponding inhomogeneous problem. The classical tool for achieving this explicitly is 
the variation-of-constants formula [30, p. 74], whose abstract formulation is given 
in Theorem 20 below. 

As usual [64], we will write &[d] for the ring of differential operators with co- 
efficients in J£\ see also Sect. 4. Let 



-n— 1 



d" 



"Co 



be a monic (i.e. having leading coefficient 1) differential operator in ,^[d] of de- 
gree n. Then we call u\,...,u n € & a fundamental system for T if it is a /T-basis 
for Ker(r), so it yields the right number of solutions for the homogeneous differen- 
tial equation Tu = 0. A fundamental system will be called regular if its associated 
Wronskian matrix 



W(ui,...,u n ) 



( «i 



u„ \ 
,1 



.(n-l) 



is invertible in ^"" x " or equivalently if its Wronskian det W(«i , . . . , u n ) is invertible 
in Of course this alone implies already that u\,...,u„ are linearly independent. 

Definition 19. A monic differential operator T e ,^[d] is called regular if it has a 
regular fundamental system. 

For such differential operators, variation of constants goes through — the canoni- 
cal initial value problem can be solved uniquely. This means in particular that regu- 
lar differential operators are always surjective. 

Theorem 20. Let , d , j) be an ordinary integro-differential algebra. Given a reg- 
ular differential operator T £ &[d] with deg T — n and a regular fundamental sys- 
tem ui,...,u n £ the canonical initial value problem 
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= em' = • 


• = em*"- 1 


> = o 



(10) 



has the unique solution 

u= 1 £u i Sd- 1 d i f (11) 

for every f G where d = &eiW(u\, . . . ,u„), anddi is the determinant of the matrix 
Wi obtained from W by replacing the i-th column by the n-th unit vector. 

Proof. We can use the usual technique of reformulating Tu = f as a system of linear 
first-order differential equations with companion matrix A G j^" x ". We extend the 
action of the operators J,<9, e componentwise to Setting now 

u = WjW" 1 / 

with / = (0,...,0,/) T G & n , we check that u G is a solution of the first- 
order system u' ~Au+f with initial condition e(m) = 0. Indeed we have u' = 
W'J W- l f + WW~ l f by the Leibniz rule and AW = W since u\ , . . . , u n are solutions 
of Tu = 0; so the differential system is verified. For checking the initial condition, 
note that eJW -1 / is already the zero vector, so we have also e(m) =0 since e is 
multiplicative. 

Writing u for the first component of m, we obtain a solution of the initial value 
problem (10), due to the construction of the companion matrix. Let us now compute 
g = W~ l f. Obviously g is the solution of the linear equation system Wg = f. Hence 
Cramer's rule, which is also applicable for matrices over rings [53, p. 513], yields 
gi as d dif and hence 

u = (Wj§) 1 = f,u i Sd- l d i f 

i=i 

since the first row of W is (mi , . . . , u n ). 

For proving uniqueness, it suffices to show that the homogeneous initial value 
problem only has the trivial solution. So assume u solves (10) with / = and choose 
coefficients ci , . . . , c„ € K such that 

m = c\ii\ H Yc n u n . 

Then the initial conditions yield e(Wc) = with c = (c\ , . . . ,c„) T G K". But we 
have also e(Wc) = (eW)c because e is linear, and detEW = E(detW) because it 
is moreover multiplicative. Since detW G ^ is invertible, eW G K" xn i s regular, so 
c=(EW) _1 = 0andM = 0. □ 
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4 Integro-Differential Operators 

With integro-differential algebras, we have algebraized the functions to be used in 
differential equations and boundary problems, but we must also algebraize the oper- 
ators inherent in both — the differential operators on the left-hand side of the former, 
and the integral operators constituting the solution of the latter. As the name sug- 
gests, the integro-differential operators provide a data structure that contains both of 
these operator species. In addition, it has as a third species the boundary operators 
needed for describing (global as well as local) boundary conditions of any given 
boundary problem for a LODE. 

Definition. The basic idea is similar to the construction of the algebra of differential 
operators i^[<3] for a given differential algebra (j^",<9). But we are now starting 
from an ordinary integro-differential algebra (^,d, J), and the resulting algebra of 
integro-differential operators will accordingly be denoted by &[d, J]. Recall that 
J£"[<9] can be seen as the quotient of the free algebra generated by d and / G 
modulo the ideal generated by the Leibniz rule df = fd +/'. For J^"[<9, J], we do 
the same but with more generators and more relations. In the following, all integro- 
differential algebras are assumed to be ordinary. 

Apart from J, we will also allow a collection of "point evaluations" as new gen- 
erators since they are needed for the specification of boundary problems. For exam- 
ple, the local boundary condition u(l) =0ona function h£^ = C°°[0, 1] gives rise 
to the functional Ei G defined by u \-¥ w(l). As one sees immediately, Ej is a 
character on J?, meaning Ei(wv) = ei(m) Ei(v) for all u,v G & . This observation 
is the key for algebraizing "point evaluations" to an arbitrary integro-differential 
algebra where one cannot evaluate elements as in C°°[0, 1]. We will see later how 
the characters serve as the basic building blocks for general local conditions like 
3u(n) — 2m(0) or global ones like J £u(£)d£. Recall that we write J?"* for the set 
of all characters on integro-differential algebra & ', In Sect. 3 we have seen that ev- 
ery integro-differential algebra («^", d, J) contains at least one character, namely the 
evaluation e = 1 — J<9 associated with the integral. Depending on the application, 
one may add other characters. 

Definition 21. Let ( d , J) be an ordinary integro-differential algebra over a field K 
and <P C The integro-differential operators J^<p[d, J] are defined as the free 
/if-algebra generated by d, and J, the "functions" / G and the characters 
<p G $U{e}, modulo the rewrite rules in Table 1. If is understood, we write 



fg ->■ f-g 


df -» fd + d-f 


in -> a-/); -/a-/) 


cpy — > \j/ 


d(p -> 


Ifd ^f-!(d-f)-(K-f)K 


<Pf ->■ (<P-/)«P 







Table 1 Rewrite Rules for Integro-Differential Operators 
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The notation U ■ /, used in the right-hand side of some of the rules above, refers 
to the action of U G ^(d, J) on a function / G J^; in particular, / ■ g denotes the 
product of two functions & . It is an easy matter to check that the rewrite rules 
of Table 1 are fulfilled in (J^",<9, J), so we may regard • as an action of J^"[<9, J] on 
Thus every element T G J? [d , J] acts as a map T : J? — > & , 

We have given the relations as a rewrite system, but their algebraic meaning is 
also clear: If in the free algebra J^(<3, J) of Definition 21 we form the two-sided 
ideal g generated by the left-hand side minus right-hand side for each rule, then 
J] = ^(^i J)/fl- Note that there are infinitely many such rules since each 
choice of e J? and <p, l// G 4> yields a different instance (there may be just 
finitely many characters in <P but the coefficient algebra & is always infinite), so g 
is an infinitely generated ideal (it was called the "Green's ideal" in [70] in a slightly 
more special setting). Note that one gets back the rewrite system of Table 1 if one 
uses the implied set of generators and a suitable ordering (see Sect. 7). 

The reason for specifying g via a rewrite system is of course that we may use it 
for generating a canonical simplifier for J^[<3, J]. This can be seen either from the 
term rewriting or from the Grobner basis perspective: In the former case, we see 
Table 1 as a confluent and terminating rewrite system (modulo the ring axioms); in 
the latter case, as a noncommutative Grobner basis with noetherian reduction (its 
elements are of course the left-hand side minus right-hand side for each rule). While 
we cannot give a detailed account of these issues here, we will briefly outline the 
Grobner basis setting since our new proof in Sect. 7 will rely on it. 

Noncommutative Grobner Bases. As detailed in Section 2, it is necessary for our 
application to deal with infinitely generated ideals and an arbitrary set of indeter- 
minates. The following description of such a noncommutative Grobner basis setting 
is based on the somewhat dated but still highly readable Bergman paper [9]; for a 
summary see [28, §3.3]. For other approaches we refer the reader to [61, 62, 85, 86]. 

Let us first recall some notions for abstract reduction relations [4]. We consider a 
relation — > C A x A for a set A; typically — > realizes a single step in a simplification 
process like the transformation of integro-differential operators according to Table 1 . 
The transitive closure of — > is denoted by -4-, its reflexive-transitive closure by 
We call a G A irreducible if there is no oq G A with a^ao; we write Aj_ for the set 
of all irreducible elements. If a Aao with ao G A^, we call «o a normal form of a, 
denoted by la = oq in case it is unique. 

If all elements are to have a unique normal form, we have to impose two con- 
ditions: termination for banning infinite reductions and confluence reuniting forks. 
More precisely, — > is called terminating if there are no infinite chains a\ — >ai — >• . . . 
and confluent \i for all a,a\,a2 G A the fork a \ A a A a? finds a reunion a \ Aao^fl2 
for some ao G A. If — > is both terminating and confluent, it is called convergent. 

Turning to noncommutative Grobner bases theory, we focus on reduction rela- 
tions on the free algebra K(X) over a commutative ring K in an arbitrary set of 
indeterminates X; the corresponding monomials form the free monoid (X). Then a 
reduction system for K(X) is a set £ C (X) x K(X) whose elements are called rules. 
For a rule a = (W,f) and monomials A,B G (X), the Zf-module endomorphism of 
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K(X) that fixes all elements of (X) except sending AWB to AfB is denoted by ^Ob 
and called a reduction. It is said to act trivially on a £ K(X) if the coefficient of AWB 
in a is zero. 

Every reduction system E induces the relation — >z C K(X) x K(X) defined by 
setting a^rzb if and only if r(a) = b for some reduction acting nontrivially on a. 
We call its reflexive-transitive closure the reduction relation induced by E, and 
we say that a reduces to b when a b. Accordingly we can speak of irreducible 
elements, normal forms, termination and confluence of E. 

For ensuring termination, one can impose a noetherian monoid ordering on (X), 
meaning a partial ordering such that 1 < A for all A £ (X) and such that B < B' 
implies ABC < AB'C for A,B,B',C £ (X). Recall that for partial (i.e. not necessarily 
total) orderings, noetherianity means that there are no infinite descending chains or 
equivalently that every nonempty set has a minimal element [8, p. 156]. Note that 
in a noetherian monoid ordering (like the divisibility relation on natural numbers), 
elements are not always comparable. 

Now if one has a noetherian monoid ordering on (X), then E will be terminating 
provided it respects < in the sense that W' < W for every rule (W,f) £ E and 
every nonzero monomial W' of /. (Let us also remark that the condition 1 < A 
from above might as well be dropped, as in [9]: The given rewrite system cannot 
contain a rule 1 — > / since then W < 1 for at least one nonzero monomial W of /, so 
1 > W > WW > • • • would yield an infinite descending chain. Such rules precluded, 
it is not stringent that constants in K be comparable with the elements in X. But 
since it is nevertheless very natural and not at all restrictive, we stick to the monoid 
orderings as given above.) 

It is typically more difficult to ensure confluence of a reduction system E. Ac- 
cording to the definition, we would have to investigate all forks a\ A a A«2, which 
are usually infinite in number. The key idea for a practically useful criterion is to 
consider only certain minimal forks (called ambiguities below, following Bergman's 
terminology) and see whether their difference eventually becomes zero. This was 
first described by Buchberger in [17] for the commutative case; see also [18, 21]. 
The general intuition behind minimal forks is analyzed in [20], where Grobner bases 
are compared with Knuth-Bendix completion and Robinson's resolution principle. 

An overlap ambiguity of E is given by a quintuple (a,T,A,B,C) with Z -rules 
a = (W,/), T = (V,g) and monomials A,B,C £ (X)\{1} such that W = AB and 
V = BC. Its associated S-polynomial is defined as fC —Ag, and the ambiguity is 
called resolvable if the S-polynomial reduces to zero. (In general one may also have 
so-called inclusion ambiguities, but it turns out that one can always remove them 
without changing the resulting normal forms [9, §5.1]. Since the reduction system 
of Table 1 does not have inclusion ambiguities, we will not discuss them here.) 

For making the connection to ideal theory, we observe that every reduction 
system E gives rise to a two-sided ideal 1^ generated by all elements W — f for 
(W,f) 6l;we have already seen this connection for the special case of the integro- 
differential operators. Note that a is equivalent to a £ 1%. 

In the given setting, the task of proving convergence can then be attacked by the 
so-called Diamond Lemma for Ring Theory, presented as Theorem 1 .2 in Bergman's 
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homonymous paper [9]; see also Theorem 3.21 in [28]. It is the noncommutative 
analog of Buchberger's criterion [18] for infinitely generated ideals. (In the version 
given below, we have omitted a fourth equivalent condition that is irrelevant for our 
present purposes.) 

Theorem 22. Let E be a reduction system for K(X) and < a noetherian monoid 
ordering that respects E. Then the following conditions are equivalent: 

• All ambiguities ofE are resolvable. 

• The reduction relation Ax is convergent. 

• We have the direct decomposition K(X) = K(X)i-\-Ie as K-modules. 

When these conditions hold, the quotient algebra K(X)/Iz may be identified with 
the K-module K(X)^, having the multiplication a-b = \,ab. 

We will apply Theorem 22 in Sect. 7 for proving that Table 1 constitutes a 
Grtibner basis for the ideal g. Hence we may conclude that &[d, j] can be identified 
with the algebra &(d, J)i of normal forms, and this is what gives us an algorith- 
mic handle on the integro-differential operators. It is thus worth investigating these 
normal forms in some more detail. 

Normal Forms. We start by describing a set of generators, which will subsequently 
be narrowed to the normal forms of J^$[<3, J"]. The key observation is that in any 
monomial we never need more than one integration while all the derivatives can be 
collected at its end. 

Lemma 23. Every integro-differential operator in ,^<p[d, J] can be reduced to a 
linear combination of monomials f(pjg\jfd', where i > and each of f,(p, J ,g, y 
may also be absent. 

Proof. Call a monomial consisting only of functions and functionals "algebraic". 
Using the left column of Table 1, it is immediately clear that all such monomials 
can be reduced to / or (p or ftp. Now let w be an arbitrary monomial in the gen- 
erators of J^&[d, J"]. By using the middle column of Table 1, we may assume that 
all occurrences of d are moved to the right, so that all monomials have the form 
w = w\ ■■■Wnd* with / > and each of W\,... ,w n either a function, a functional 
or J. We may further assume that there is at most one occurrence of J among the 
w\,... ,w„. Otherwise the monomials w\ ■ ■ ■ w n contain JwJ, where each w = ftp is 
an algebraic monomial. But then we can reduce 

M=a/<p)j=a-/)<pj 

by using the corresponding rule of Table 1 . Applying these rules repeatedly, we 
arrive at algebraic monomials left and right of J (or just a single algebraic monomial 
if J is absent). □ 

In TH30REMV, the integro-differential operators over an integro-differential al- 
gebra & of coefficient functions are built up by FreelntDif f Op[J?,K]. This func- 
tor constructs an instance of the monoid algebra with the word monoid over the 
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infinite alphabet consisting of the letters d and J along with a basis of & and with 
all multiplicative characters induced by evaluations at points in K: 



Definition ["IntDiff Op", any[:r, K] , 
IntDiffOp[r, K] = where[« = FreelntDif fOp[r, K] , § - GreenSystem[r, K] 
QuotAlg[GBNF[«, g]]] 



In this code fragment, the GreenSystem functor contains the encoding of the afore- 
mentioned rewrite system (Table 1), here understood as a noncommutative Grobner 
basis. Normal forms for total reduction modulo infinite Grobner bases are created by 
the GBNF functor, while the QuotAlg functor constructs the quotient algebra from 
the corresponding canonical simplifier (see Sect. 2 for details). For instance, multi- 
plying the integral operator J by itself triggers an application of the Baxter rule: 



Here integral operators are denoted by A, following the notation in the older imple- 
mentation [70]. 

We turn now to the normal forms of boundary conditions. Since they are intended 
to induce mappings & — > K, it is natural to define them as those integro-differential 
operators that "end" in a character (p G <P. For example, if <p is the point evaluation 
Ei considered before, the composition ei<3 describes the local condition u'(l) = 
0, the composition Ei j the global condition jlu(£)d£ = 0. In general, boundary 
conditions may be arbitrary linear combinations of such composites; they are known 
as "Stieltjes conditions" in the literature [15, 16]. 

Definition 24. The elements of the right ideal 



are called boundary conditions over & , 

It turns out that their normal forms are exactly the linear combinations of local 
and global conditions as in the example mentioned above. As a typical representative 
over & = C°°(R), one may think of an element like 



Compute [((l, ("/")))* {{l, ("J")))] // FormatP 



-A X + X A 



<P) = <P.J? iP [dJ] 



Eo<3 2 + 3 Ejr — 2 E2jt J sinx, 
written as u"(0) + 3u(n) — sin B,u(t,)dB, in traditional notation. 



Proposition 25. Every boundary condition of\<P) has the normal form 




(12) 



with only finitely many a^i £ K and /<p £ & nonzero. 
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Proof. By Lemma 23, every boundary condition of \<P) is a linear combination of 
monomials having the form 

w = Xf<Phvd i or w = xf<pd' (13) 

where each of /, g, (p , yi may also be missing. Using the left column of Table 1 , the 
prefix xf*P can be reduced to a scalar multiple of a functional, so we may as well 
assume that / and <p are not present; this finishes the right-hand case of (13). For 
the remaining case w = %§gyd\ assume first that y/ is present. Then we have 

x(Sgv)=x(S-g)Y= (xS-8)xy = (xj-g) v> 

so w is again a scalar multiple of yfd', and we are done. Finally, assume we have 
w = X§gd'- ^ i = 0' m is is already a normal form. Otherwise we obtain 

w = X (Jgd) d i - 1 = (x-g)xd i - 1 -xSg'd i - 1 -(*-g)*d l -\ 

where the first and the last summand are in the required normal form, while the 
middle summand is to be reduced recursively, eventually leading to a middle term 
in normal form ±xjg'd° = ±xjg'. □ 

Most expositions of boundary problems — both analytic and numeric ones — 
restrict their attention to local conditions, even more specifically to those with just 
two point evaluation (so-called two-point boundary problems). While this is doubt- 
less the most important case, there are at least three reasons for considering Stieltjes 
boundary conditions of the form (12). 

• They arise in certain applications (e.g. heat radiated through a boundary) and in 
treating ill-posed problems by generalized Green's functions [70, p. 191]. 

• As we shall see (Sect. 5), they are needed for factoring boundary problems. 

• Their algebraic description as a right ideal is very natural. 

Hence we shall always mean all of \<P) when we speak of boundary conditions. 

Let us now turn to the other two ingredients of integro-differential operators: We 
have already mentioned the differential operators J£"[<9], but we can now see them as 
a subalgebra of JP<p[d, J"]. They have the usual normal forms since the Leibniz rule 
is part of the rewrite system. Analogously, we introduce the subalgebra of integral 
operators generated by the functions and J. Using Lemma 23, it is clear that the 
normal forms of integral operators are & itself and linear combinations of fjg, and 
the only rule applicable to them is J/ J —>•••■ in Table 1 . Since we have already 
included & in J?[d], we introduce as the i^-bimodule generated by J so that 
it contains only the monomials of the form fjg. 

Finally, we must consider the two-sided ideal (<P) of &®\d, J] generated by <P 
whose elements are called boundary operators. A more economical description of 
(4>) is as the left J^-submodule generated by |<f>) because by Lemma 23 any wxw 
with w, w G & [d, j] can be reduced to f(pjgyrd'xw. Note that (<P) includes all finite 
dimensional projectors P along Stieltjes boundary conditions. Any such projector 
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can be described in the following way: If ui, . . . ,u n G & and j3i , . . . ,/3„ G \<P) are 
biorthogonal (meaning Pi{uj) = 5;/), then 

P = X>/3,:^^^ (14) 

is the projector onto [«i , . . . , u n ] whose kernel is the subspace of all such that 

/3(w) =0 for all /3 G [j3i , . . . , j3„] . See for example [50, p. 71] or [66] for details. Note 
that all elements of (<P) have the normal form (14), except that the (uj) need not be 
biorthogonal to the (/?;). 

We can now characterize the normal forms of J^>[<3, J] in a very straightforward 
and intuitive manner: Every monomial is either a differential operator or an integral 
operator or a boundary operator. Hence every element of ^<$[d, j] can be written 
uniquely as a sum T + G + B, with a differential operator T G J^[<3], an integral 
operator G G .i^[f], and a boundary operators G (<£). 

Proposition 26. For an integro-differential algebra a«c/ characters <P C j?*, we 
/zave f/ie direct decomposition J^>[<9, J] = J^"[<9] + + 

Proof. Inspection of Table 1 confirms that all integro-differential operators having 
the described sum representation T + G + P are indeed in normal form. Let us now 
prove that every integro-differential operator of ,^<p [<3 , J] has such a representation. 
It is sufficient to consider its monomials w. If w starts with a functional, we ob- 
tain a boundary condition by Proposition 25; so assume this is not the case. From 
Lemma 23 we know that 

w = f(pj'g\j/d' or w = f(pd', 

where each of (p,g, y/ may be absent. But w G (<P) unless <p is absent, so we may 
actually assume 

w = fjgyd' or w=fd*. 

The right-hand case yields w G i^[<3]. If is present in the other case, we may 
reduce jgy/ to (J ■ g) \]/, and we obtain again w G (<£>). Hence we are left with w = 
fjgd', and we may assume ; > since otherwise we have w G i^[f] immediately. 
But then we can reduce 

W = f(Sgd)d i - l =f(g-S(d-g)-(s-g)z)d i - 1 

= (fg)d i - l -fS(d-g)d- 1 -(K-g)fKd i -\ 

where the first term is obviously in ,^[d] and the last one in (<£). The middle term 
may be reduced recursively until the exponent of d has dropped to zero, leading to 
a term in &\ J]. □ 

We can observe the direct decomposition J£<j>[<9, J] = ^[d] + J^"[f] + (4>) in the 
following sample multiplication of J<9 and ddxe x J: 
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Compute[((l, ("J", "B")))*{{1, {"3", "3", <»|"l", <1, 1», "J")))] // FormatP 

-2E + 2e x + 2e x x + 2e x A + e x xA + e x xD 

As in the previous computation, A stands for the integral J, moreover D for the 
derivation d, and E for the evaluation. As we can see, the sum is composed of one 
differential operator (the last summand), two integral operators (in the middle), and 
three boundary operators (the first summands). Observe also that the input operators 
are not in normal form but the output operator is. 

Basis Expansion. Regarding the canonical forms for <^{d, J"], there is one more 
issue that we have so far swept under the rug. The problem is that in the current 
setup elements like xd + 3x 2 <9 and (x + 3x 2 )d are considered distinct normal forms. 
More generally, if / + g = h holds in & , there is no rule that allows us to recognize 
that / + g G &\d, j] and h G J£"[<9, J] are the same. Analogously, if Xf = g holds 
in & with X G K, then Xf and g are still considered to be different in J^"[<3, J]. A 
slightly less trivial example is when / = (cosjc)(cos 2 x 2 ) and g = — (sinx)(sinx 2 ) so 
that h = cos (x + x 2 ). What is needed in general is obviously a choice of basis for 
J£\ But since such a choice is always to some degree arbitrary, we would like to 
postpone it as much as possible. 

An unbiased way of introducing all /f-linear relations in & is simply to collect 
them in all in the two-sided ideal 

l=(f+g-h,Xf-g\f+g = hmdXf = gm &)C&(dJ), 

which we shall call the linear ideal. Since [ + g corresponds to a unique ideal [ in 
& [d , J] , the necessary refinement of & [d , J] can now be defined as 

,? # ldJ}=,ndJ}/l = ^{dJ)/(l + Q) 

whose elements shall be called expanded integro-differential operators. Note that [ 
is really the "same" ideal as [ except that now f,g,h,f,g G ^[d, J"]. By the isomor- 
phism above, coming from the Third Isomorphism Theorem [31, Thm 1.22], we can 
think of J^ # [<9, J] in two ways: Either we impose the linear relations on J^[d, J] or 
we merge them in with the Green's ideal — let us call these the a-posteriori and the 
combined approach, respectively. 

For computational purposes, we need a ground simplifier on the free algebra [70, 
p. 525], which we define here as a Zf-linear canonical simplifier for J^(<9,J)/[. 
Since all reduction rules of Table 1 are (bi)linear in /, g G any ground simplifier 
descends to a canonical simplifier a on [<?,/]. In our implementations, a always 
operates by basis expansion (see below), but other strategies are conceivable. We 
can apply a either a-posteriori or combined: 

• In the first case we apply ffasa postprocessing step after computing the normal 
forms with respect to Table 1. We have chosen this approach in the upcoming 
Maple implementation [49]. 

• In the combined approach, a may be used at any point during a reduction along 
the rules of Table 1 . It may be more efficient, however, to use a on the rules 
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themselves to create a new reduction system on J^"(<?, J); see below for an ex- 
ample. We have taken this approach in our earlier implementation [69, 70] and 
in the current implementation. 

Generally the first approach seems to be superior, at least when a tends to create 
large expressions that are not needed for the rewriting steps of Table 1 ; this is what 
typically happens if the ground simplifier operates by basis expansion. 

Assume now we choose a -basis (bj)iei of & . If (oj)/e? is the dual basis, we 
can describe the linear ideal more economically as 

[=(/-!>(/)&, | /e«^), 

16/ 

so the linear basis (bi)i e i gives rise to an ideal basis for L Its generators /— £ ; • • • can 
be oriented to create a ground simplifier a : / i-> £ ; ■ • • effecting basis expansion. 

If one applies now such a ground simplifier coming from a basis (£>,), G / in the 
combined approach, one can restrict the generators of &(d, j) to basis elements 
bi G & rather than all / G and the reduction rules can be adapted accordingly. 
For example, when & contains the exponential polynomials, the Leibniz rule df —> 
fd + (d ■ f) gets instantiated for / = xe x as d{xe x ) — > (xe x )d + e x + x, where the 
right-hand side now has three instead of two monomials! This is why the choice of 
basis was built into the definition of the precursor of J^" # [<3 , J] as in [70] . 

Before leaving this section on integro-differential operators, let us mention some 
interesting related work on these objects, carried out from a more algebraic view- 
point. In his papers [5, 6], Bavula establishes an impressive list of various (notably 
ring-theoretic) properties for algebras of integro-differential operators. The setup 
considered in these papers is, on the one hand, in many respects more general since 
it deals with partial rather than ordinary differential operators but, on the other hand, 
the coefficients are restricted to polynomials. 

Seen from the more applied viewpoint taken here, the most significant difference 
is the lack of multiple point evaluations (and thus boundary conditions). Apart from 
these obvious differences, there is also a somewhat subtle difference in the meaning 
of e = 1 — Jo<3 that we have tried to elucidate in a skew-polynomial setting [67]. 
The upshot is that while our approach views e as a specific evaluation (the proto- 
typical example is given after Definition 2), it does not have a canonical action in 
V. Bavula's setting (and neither in our skew-polynomial approach). This is a subtle 
but far-reaching difference that deserves future investigation. 



5 Applications for Boundary Problems 

In this section we combine the tools developed in the previous sections to build an 
algorithm for solving linear boundary problems over an ordinary integro-differential 
algebra; see also [73] for further details. We also outline a factorization method 
along a given factorization of the defining differential operator applicable to bound- 
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ary problems for both linear ordinary and partial differential equations; see [66] in 
an abstract linear algebra setting and [74] for an overview. 

For motivating our algebraic setting of boundary problems, let us consider our 
standard example of an integro-differential algebra (j^", d, J) with the integral oper- 
ator 

}■■/■ ■ 

for & = C°°[0, 1]. The simplest two-point boundary problem reads then as follows: 
Given / 6 find such that 



» — j i 

«(0) = k(1) = 0. 



In this and the subsequent examples, we let D and A denote respectively the deriva- 
tion d and the integral operator J. Moreover, we denote by L the corresponding 
evaluation e, which is the character 

L:/->/(0). 

To express boundary problems we need additionally the evaluation at the endpoint 
of the interval 

R:f^f(l). 

Note that u is annihilated by any linear combination of these functionals so that 
problem (15) can be described by the pair (D 2 , [L,R]), where [L,R] is the subspace 
generated by L, R in the dual space J^"* . 

The solution algorithm presupposes a constructive fundamental system for the 
underlying homogeneous equation but imposes no other conditions (in the liter- 
ature one often restricts to self-adjoint and/or second-order boundary problems). 
This is always possible (relative to root computations) in the important special case 
of LODEs with constant coefficient. 

The Solution Algorithm. In the following, we introduce the notion of boundary 
problem in the context of ordinary integro-differential algebras. Unless specified 
otherwise, all integro-differential algebras in this section are assumed to be ordinary 
and over a fixed field K. 

Definition 27. Let (J^",<9,J) be an ordinary integro-differential algebra. Then a 
boundary problem of order n is a pair (T,£§), where T G ^[d] is a regular differ- 
ential operator of order n and 33 C | J?*) is an n-dimensional subspace of boundary 
conditions. 

Thus a boundary problem is specified by a differential operator T and a bound- 
ary space 33 = [j3i, . . . ,J3„] generated by linearly independent boundary conditions 
J3i , . . . , J3„ £ | & % ) . In traditional notation, the boundary problem (T,3§)is then given 
by 
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Tu = f, 

j3\u = ■■■ = j8„« = 0. 



(16) 



For a given boundary problem, we can restrict to a finite subset <P C with the 
consequence that all subsequent calculations can be carried out in &&[d, J] instead 
of &[d, J] . We will disregard this issue here for keeping the notation simpler. 

Definition 28. A boundary problem (T, 33) is called regular if for each / G & there 
exists a unique solution u £ & in the sense of (16). 

The condition requiring T to have the same order as the dimension of S3 in Def- 
inition 27 is only necessary but not sufficient for ensuring regularity: the boundary 
conditions might collapse on Ker(T). A simple example of such a singular bound- 
ary problem is (— D 2 , [LD 7 RD]) using the notation from before; see also [70, p. 191] 
for more details on this particular boundary problem. 

For an algorithmic test of regularity, one may also apply the usual regularity 
criterion for two-point boundary problems, as described in [66]. Taking any funda- 
mental system of solutions u\,...,u n for the homogeneous equation, one can show 
that a boundary problem (T, S3) is regular if and only if the evaluation matrix 



/j3i(«i) ■•■ /3i("«) N 
£(«)= : ■-. : 
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is regular. 

For a regular boundary problem (T, 33), we can now define its Green 's operator 
G as the linear operator mapping a given forcing function / G & to the unique 
solution u G & of (16). It is characterized by the identities 

TG=\ and ha(G)=^ ± , 

where B3 X = {u G & \ J3(w) = for all j3 G £§} is the subspace of all "functions" 
satisfying the boundary conditions. We also write 



G={T,BS) 



-l 



for the Green's operator of (T, BS). 

The investigation of singular boundary problems (i.e. non-regular ones) is very 
enlightening but leads us too far afield; we shall investigate it at another junction. 
Let us just mention that it involves so-called modified Green's operators and func- 
tions [80, p. 216] and that is paves the way to an interesting non-commutative analog 
of the classical Mikusiriski calculus [60]. 

We will now recast Theorem 20 in the language of Green's operators of ini- 
tial value problems. Given a regular differential operator T of order n, the theorem 
implies that the initial value problem (T, [e, e<9, . . . , e<9" -1 ]) is regular. We call its 
Green's operator the fundamental right inverse of T and denote it by T*. 
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Corollary 29. Let («^*, d, j) be an ordinary integro-differential algebra and let T G 
&[d] be a regular differential operator of order n with regular fundamental system 
Ml, . . . ,U n . Then its fundamental right inverse is given by 

n 

T* = Y,Ui$d- 1 d i G &[dj], (17) 

where d,d\,...,d n are as in Theorem 20. 

Before turning to the solution algorithm for boundary problems, let us also men- 
tion the following practical formula for specializing Corollary 29 to the important 
special case of LODEs with constant coefficients, which could also be proved di- 
rectly e.g. via the Lagrange interpolation formula. For simplicity, we restrict our- 
selves to the case where the characteristic polynomial is separable. 

Corollary 30. Let (J?,<3, J) be an ordinary integro-differential algebra and con- 
sider the differential operator T = (<9 — Aj) • • • (<3 — A„) G &[d] with Ai,. . . ,A„ GAT 
mutually distinct. Assume each u' = A,- u, e • u = 1 has a solution u = e^' x G & with 
reciprocal m _i = e~^' x G Then we have 

i=i 

where yL[ = (A,- — Ai) • • • (A; - A,--i)(A; - Aj+i) • • • (A,- - A„). 

Proof. Let us write V for the n x n Vandermonde determinant in Ai , . . . , A„ and 
Vi for the (n — 1) x (n — 1) Vandermonde determinant in Ai , . . . , A,_ i , A, + 1 , . . . , A„ . 
Evaluating the quantities of (17), one sees immediately that 

d=e (M+-l„)xy an( j di= ^_ l y+i e (l l +...+X i _ l +l i+l +-+X n )x Vh 

Hence we have dj/d = (— l) n+l e~^ x Vj/V. Using the well-known formula for the 
Vandermonde determinant, one obtains d\jd = \iie~^> x , and now the result follows 
from Corollary 29. □ 

Summarizing our earlier results, we can now give a solution algorithm for com- 
puting G = (T,38)~ l , provided we have a regular fundamental system ui, . . . ,u„ for 
Tu = and a /iT-basis J3i, . . . , J3„ for 38. The algorithm proceeds in three steps: 

1 . Construct the fundamental right inverse r * G J? [<3 , J ] as in Corollary 29 . 

2. Determine the projector P = Y!i=\ u ifii £ ^[d, J] as in (14). 

3. Compute G = (l-P)T* e &[d,J]. 

Theorem 31. The above algorithm computes the Green 's operator G G ^[d, J] for 
any regular boundary problem (T,3§). 



Proof. See [73]. □ 
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The computation of Green 's operators for boundary problems for ODEs using 
the above algorithm takes on the following concrete form in TH30REMV code. 

GreensOp[F, 2] = [l - Proj [2, F] ] * RiqhtInv[F] 

p I* * p I * p 

Here SB is a basis for the boundary space and F a regular fundamental system. 
Let us consider again example (15): Given / £ & ', find such that 



«(0) = k(1) = 0. 



The Green's operator G of the boundary problem can be obtained by our implemen- 
tation via the following computation 

G = GreensOp[D 2 , <L, R>] 

-A x-xB + xAx + xBx 



where we use the notation from before: An = Jg u(%) d% , Lu = u(Q), Ru = u(l) and 
in addition, Bit = u(£)d£. The corresponding Green's function is computed in 
an immediate postprocessing step: 

GreensFct [G] 

. „ r-f + xf «= e<x 
g[x, a = 

L-x + xi; <= X < .* 

As noted in [70], the Green's function provides a canonical form for the Green's 
operator. Moreover, one can obtain the function u(x) and thus solve the boundary 
problem through knowledge of the Green's function in the following identity: 

u{x) = Gf{x)= f ' g(x,£,)f{S)d$. 
Jo 

By replacing the Green's function obtained above in the latter integral we obtain 

u (x) = (x - 1 ) r d^ +x C (§ - 1 )M ) dt; . 

Jo Jx 

Furthermore, we can look at some specific instances of the forcing function f(x). 
Let us first consider the simple example f(x) = x. By an immediate calculation, we 
obtain the expression for the action of G on f(x), which is u(x): 

GreensOpAct [G, x] 

- + ^: +X 3 f-ii 

6 3 I, 6 J 

The expression for the solution function u(x) can easily become more complicated, 
as we can see in the next example, where we consider the instance 



f(x) =e 2A " + 3x 2 sinx 3 . 
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Relying on Mathematica for handling symbolic integration, we obtain: 



GreensOpAct[G, e 2 x + 3 x 2 Sin [x] 3 ] 



- 1 e 2 x x € 
4 4 4 


2 x f-31 

— + — e 2 * x + e 2 * x 2 

4 { 2 ) 


1 

+ ( - 9) x Cos [1] + - xCos[3] + 
9 


-18) xCos[x] + 27x 2 Cos[x 


9 1 
- x 3 Cos [x] + 

2 i 


-91 2 

x 4 Cosfx] + — x Cos 

2 j 9 


[3 x] + |—^J x 2 Cos [3 x] + I— -J 


x 3 Cos [3 x] + 


1 

- x* Cos [3 x] + 
2 


45 1 

— xSinfl] + — xSinf 
4 36 


27 Sin [x] 

3] + + (-27) xSinfx 

2 


* x 2 Sin[x] ♦ 


27 

— x 3 Sin [x] + 

2 


(-It 1 

— Sin[3 x] + — x Sin 
Us j 9 


5 / -1 \ 

[3 x] + — x 2 Sin [3 x] + \ x' 

12 [2 


Sin[3 x] 



As a last example, let us consider f(x) = sinsinx. As we can notice below, it cannot 
be integrated with Mathematica: 

GreensOpAct[G, Sin[Sin[x]]] 

J (-x.xE) Sin[Sin[f]] + (-ftx£) Sin[Sin[f]] d£ 

In order to carry out the integrals involved in the application of the Green's oper- 
ator to a forcing function, one can use any numerical quadrature method (as also 
available in many computer algebra systems). 



Composing and Factoring Boundary Problems. In the following, we discuss the 
composition of boundary problems corresponding to their Green's operators. We 
also describe how factorizations of a boundary problem along a given factoriza- 
tion of the defining operator can be characterized and constructed. We refer again 
to [73, 66] for further details. We assume that all operators are defined on suitable 
spaces such that the composition is well-defined. It is worth mentioning that the fol- 
lowing approach works in an abstract setting, which includes in particular boundary 
problems for linear partial differential equations (LPDEs) and systems thereof; for 
simplicity, we will restrict ourselves in the examples to the LODE setting. 

Definition 32. We define the composition of boundary problems (T\,,9§\) and (T 2 ,3§ 2 ) 
by 

(7"i , ®x) o {Jz,3h) = (TiT 2 , ®x ■ T 2 + 3h) ■ 

So the boundary conditions from the first boundary problem are "translated" by 
the operator from the second problem. The composition of boundary problems is 
associative but in general not commutative. The next proposition tells us that the 
composition of boundary problems preserves regularity. 

Proposition 33. Let (Ti,3§i) and {T2,S§2) be regular boundary problems with 
Green's operators G\ and G 2 - Then (T\,£$\) o (72, ,^2) is regular with Green's op- 
erator G 2 G\ so that 

((7i,#i)o (r 2 ,^ 2 ))- 1 = (T 2 ,mr 1 °(T h m)- 1 . 



The simplest example of composing two boundary (more specifically, initial 
value) problems for ODEs is the following. Using the notation from before, one 
sees that 
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(D, [L])o (D, [L]) = (D\ [LD] + [£]) = (D 2 , [L,LD]). 

Let now (T,3§) be a boundary problem and assume that we have a factorization 
T = T\T2 of the defining differential operator. We refer to [73, 66] for a characteri- 
zation and construction of all factorizations 

(T,&) = (Ti,&i) (T 2t &2) 

into boundary problems. In particular, if (T,3&) is regular, it can be factored into 
regular boundary problems: the left factor (T\,3S{) is unique, while for the right 
factor (72, &8z) we can choose any subspace 8$i < £$ that makes the problem regu- 
lar. We can compute the uniquely determined boundary conditions for the left factor 
by SS\ = SB ■ G2, where G2 is the Green's operator for some regular right factor 
(T2,&2). By factorization, one can split a problem of higher order into subproblems 
of lower order, given a factorization of the defining operator. For algorithms and 
results about factoring ordinary differential operators we refer to [64, 78, 84]. 

Given a fundamental system of the differential operator T and a right inverse of 
T2, one can factor boundary problems in an algorithmic way as shown in [66] and 
in an integro-differential algebra [73]. As described in [74], we can also compute 



boundary conditions , 



<, 



' such that ( T2 , S&2 ) is a regular right factor, given only a 



fundamental system of T2. The unique left factor can be then computed as explained 
above. This allows us to factor a boundary problem if we can factor the differential 
operator and compute a fundamental system of only one factor. The remaining lower 
order problems can then also be solved by numerical methods. 

Here is how we can compute the boundary conditions of the left and right fac- 
tor problems for the boundary problem (D 2 ,[L,R)) from previous example (15), 
along the trivial factorization with T\ = T2 = D. The indefinite integral A = Jq is the 
Green's operator for the regular right factor (D, [L] ). 

Fact[D, D, <L, R>, <R>] 

<<A + B>, (L>> 

This factorization reads in traditional notation as 



u'=f 
■1 



u' = f 
tf(0) = 



u"=f 

«(0) =m(1) =0 



Note that the boundary condition for the unique left factor is an integral (Stieltjes) 
boundary condition. 

We consider as a second example the fourth order boundary problem [73, Ex. 33]: 



u"" + 4u = f, 

m(0) =u(\) =m'(0) =m'(1) =0. 



(18) 



Factoring the boundary problem along D 4 + 4 = (D 2 — 2i) (D 2 + 2i), we obtain the 
following boundary conditions for the factor problems. 
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Fact[D 2 -2i, D 2 + 2i, <L, R, LD, R D) , (a'" 1 * 1 " 1 , s 11 "* 1 " 

( ( A e (Complex! -1,1] I » + B e (Complex! -1,1] ) ^ fl e (Complex [ 1, -1 ] | ; 




6 Integro-Differential Polynomials 

In this section, we describe the algebra of integro-differential polynomials [72] ob- 
tained by adjoining an indeterminate function to a given integro-differential algebra 
(J^",<9, J). Intuitively, these are all terms that one can build using an indeterminate 
u, coefficient functions / G & and the operations +,-,d, J, identifying two terms if 
one can be derived from the other by the axioms of integro-differential algebras and 
the operations in J? ', A typical term for (K[x] , <3, J) looks like this: 



From the computational point of view, a fundamental problem is to find a canon- 
ical simplifier (see Sect. 2) on these objects. For example, the above term can be 
transformed to 



by the Baxter axiom and the /if-linearity of the integral. 

As outlined in the next subsection, a notion of polynomial can be constructed for 
any variety in the sense of universal algebra. (In this general sense, an algebra is 
a set with an arbitrary number of operations, and a variety is a collection of such 
algebras satisfying a fixed set of identities. Typical examples are the varieties of 
groups, rings, and lattices.) 

For sample computations in the algebra of integro-differential polynomials, we 
use a prototype implementation of integro-differential polynomials, based on the 
TH30REMV functor mechanism (see Sect. 2). 

Polynomials in Universal Algebra. In this subsection, we describe the idea of the 
general construction of polynomials in universal algebra [45]. We refer to [54] for 
a comprehensive treatment; see also the surveys [1, 27]. For the basic notions in 
universal algebra used below, see for example [4] or [33, Ch. 1]. 

Let V be a variety defined by a set S of identities over a signature E. Let A be a 
fixed "coefficient domain" from the variety "V , and let X be a set of indeterminates 
(also called "variables"). Then all terms in the signature E with constants (hence- 
forth referred to as coefficients) in A and indeterminates in X represent the same 
polynomial if their equality can be derived in finitely many steps from the identities 
in <§ and the operations in A. The set of all such terms S?z (A UX) modulo this con- 
gruence = is an algebra in "f , called the polynomial algebra for "f in X over A and 
denoted by A^[X]. 



(4mm' J> + 3)m' 3 )(m' Jm" 2 ) + JVV' 5 J> 2 +5x)u 3 u' 2 ju 



4mm' 2 Jxm' 3 Jm" 2 +4mm ,2 Jm" 2 Jxm' 3 + 12mm' 2 Jm' 3 Jm" 2 + 12mm' 2 Jm" 2 Jm' 

+ jx 6 uu" 5 $x 2 u\' 2 ju + 5$x 6 uu" 5 $i 



xm 3 m' 2 Jm. 
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The polynomial algebra Ay [X] contains A as a subalgebra, and A U X is a gener- 
ating set. As in the case of polynomials for commutative rings, we have the substitu- 
tion homomorphism in general polynomial algebras. Let B be an algebra in "¥ , Then 
given a homomorphism q>\ : A — > B and a map q>2 : X — >• B, there exists a unique 
homomorphism 

<p: A r [X] ->B 

such that (p(a) = (pi (a) for all a E A and <p(jc) = <P2(*) f° r a ll x E X. So in a categor- 
ical sense the polynomial algebra Ay [X] is a free product of the coefficient algebra 
A and the free algebra over X in see also [1]. 

For computing with polynomials, we will construct a canonical simplifier on 
Ay [X] with associated system of canonical forms M. As explained before (Sect. 2), 
the canonical simplifier provides for every polynomial in Ay[X], represented by 
some term T, a canonical form R E & that represents the same polynomial, with 
different terms in representing different polynomials; see also [54, p. 23]. 

The set ffl must be large enough to generate all of Ay[X] but small enough to 
ensure unique representatives. The latter requirement can be ensured by endowing 
a given set of terms with the structure of an algebra in the underlying variety. 

Proposition 34. Let "V be a variety over a signature E, let A be an algebra in "V 
and X a set of indeterminates. If 8% C &z{A UX) is a set of terms with AUX C 3% 
that can be endowed with the structure of an algebra in "f ', then different terms in 
8% represent different polynomials in Ay [X], 

Proof. Since 8i can be endowed with the structure of an algebra in the variety "V 
and A UlCil, there exists a unique substitution homomorphism 

<p:Ay\X]^M 

such that (p(a) =a for all a E A and <p(x) = x for all x E X. Let 

k: 8f^A r [X] 

denote the restriction of the canonical map associated with =. Then we have (p o 
Tt{R) = R for all R E 8%, so n is injective, and different terms in ffl indeed represent 
different polynomials. □ 

As a well-known example, take the polynomial ring R[x] in one indeterminate x 
over a commutative ring R, which is A y [X] for A = R and X = {x} with "f being the 

variety of commutative unital rings. The set of all terms of the form a n x n H h «o 

with coefficients a; E R and a„ ^ together with is a system of canonical forms 
for R [x] . One usually defines the polynomial ring directly in terms of these canonical 
forms. Polynomials for groups, bounded lattices and Boolean algebras are discussed 
in [54] along with systems of canonical forms. 

Differential Polynomials. For illustrating the general construction described above, 
consider the algebra of differential polynomials over a commutative differential K- 
algebra (J?,<3) in one indeterminate function u, usually denoted by ,^{u}. Clearly 
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this is Ay[X] for A = & and X = {u} with f being the variety of differential K- 
algebras. Terms are thus built up with the indeterminate u, coefficients from & and 
the operations +,-,d; a typical example being 

d 2 (f lU 2 + U )d(f 2U 3 ) + d 3 (f }U ). 

By applying the Leibniz rule and the linearity of the derivation, it is clear that every 
polynomial is congruent to a A'-linear combination of terms of the form 

fjj«i, (19) 

i=0 

where / S the notation u„ is short for d n (u), and only finitely many /3, € N are 
nonzero. In the following, we use the multi-index notation fur for terms of this 
form. For instance, 

M (l,0,3,2) ig 

the multi-index notation for u(u") 3 (u"') 2 . The order 
of a differential monomial is given by the highest derivative appearing in ift or 
-ooif J3 =0. 

Writing & for the set of all K-Ymzax combinations of terms of the form (19), we 
already know that every polynomial is congruent to a term in M. When & = K[x], 
a typical element of & is given by 

(3x 3 + 5x) + 7*V 2 - ai ) + 2xu^\ 

To show that 3% is a system of canonical forms for ^{u}, by Proposition 34 it 
suffices to endow ffl with the structure of a commutative differential algebra. As 
a commutative algebra, ffl is just the polynomial algebra in the infinite set of in- 

determinates koj m 1) m 2) For defining a derivation in a commutative algebra, by 

the Leibniz rule and /^-linearity, it suffices to specify it on the generators. Thus 
ffl becomes a differential algebra by setting = itk+i- One usually defines the 

differential polynomials directly in terms of these canonical forms, see for exam- 
ple [48]. 

Integro Polynomials. We outline the integro polynomials over a Rota-Baxter algebra 
as in Definition 8. This is related to the construction of free objects in general Rota- 
Baxter algebras; we refer to [41] for details and references. By iterating the Baxter 
axiom (8), one obtains a generalization that is called the shuffle identity on &\ 

( J/i J • • • J/») • (Jsi J • • • Jg») = £ JW • • • JW« (20) 

Here the sum ranges over all shuffles of (/i, .. . ,/,„) and (gi,... ,g n )l see [65, 68, 76] 
for details. The sum consists of (" ! ^") shuffles, obtained by "shuffling" together 
J/i J ' ' ' J/"' an d Ssi J " ' §8n as words over the letters jf and jgj, such that the 
inner order in the words is preserved. For instance, we have 



(J/1J/2) • (/si) = J/i Sfitei + J/i hijfi + JgiSfiSfi- 
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for the simple m = 2,n = 1 case. 

The integro polynomials over are defined as Ay[X] for A = j£" and X = {«} 
with "f being the variety of Rota-Baxter algebras over K. The full construction of 
the canonical forms for integro polynomials is included in the following subsection. 
But it is clear that by expanding products of integrals by the shuffle identity, every 
integro polynomial is congruent to a /if-linear combination of terms of the form 

/n*j/i «**/. ••;/„,«*»■ (2i) 

with f,fi,...,f m £ & an d k,k\, . . . ,k m £ N. However, they cannot be canonical 
forms, since terms like j(f + g)u and jfu + jgu or jXfu and A jfu represent the 
same polynomials. 

Writing ffi for the set of all Zf-linear combinations of terms of the form (21), the 
multiplication of two elements of S% can now be defined via (20) as follows. Since 
the product of (21) with another term gu 1 §g\u tl / • • • §g„u l " should clearly be given 
by fgu k+, (jfiu kl j---jf m u km )(jgu'giu' 1 j--- f n u'"), it remains to define the so- 
called shuffle product on integral terms (those having the form (21) with / = 1 and 
k = 0). This can be achieved immediately by using (20) with fu k ' and gju'i in place 
of f and gj, respectively. It is easy to see that the shuffle product is commutative 
and distributive with respect to addition. 

The shuffle product can also be defined recursively [68]. Let J and / range over 
integral terms (note that 1 is included as the special case of zero nested integrals). 
Then we have 

(//«*/) ■ (Jfu 1 J) = (jfu k ) U J ■ ( tfjj) + (]fu k ) U ( $fu k J) ■ J, (22) 

where U : 3% x — > denotes the operation of nesting integrals (with the under- 
standing that ■ binds stronger than U), defined on basis vectors by 

J>i $F m u Jd J ... JG„ = .|7,.r • • • J7i,,|V; : J. • • JG„, (23) 

and extended bilinearly to all of 3?. Here F, and Gj stand for fur' and gjU 1 ' , respec- 
tively. For example, $F\ J7*2 an d jGi can be multiplied as 

(pi) U (JF 2 ) ■ (/<?!)+ (/d) U 1 • (pijF 2 ) = (/Fx) U (/F 2 /d +/G!/F 2 ) 
+ CfGi) U (/Fi/F 2 ) = /Fx/Fi/Gj +/F 1 /Gi/F 2 + /G 1 /F 1 /F 2j 

analogous to the previous computation. 

Representing Integro-Differential Polynomials. In the following, we describe in 
detail the universal algebra construction of the integro-differential polynomials and 
their canonical forms. We refer to [39, 40] for the related problem of constructing 
free objects in differential Rota-Baxter algebras. We consider the variety of integro- 
differential algebras. Its signature L contains: the ring operations, the derivation d, 
the integral /, the family of unary "scalar multiplications" {•\)xek> an d f° r conve- 
nience we also include the evaluation e. The identities S are those of a Zf-algebra, 
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then ^-linearity of the three operators d, J, e, the Leibniz rule (2), the section axiom 
(3), the Definition 2 of the evaluation, and the differential Baxter axiom (6). 

Definition 35. Let (jF, d , J) be an integro-differential algebra. Then the algebra of 
integro-differential polynomials in u over denoted by ,^{u} in analogy to the 
differential polynomials, is the polynomial algebra Ay [AT] for A = & and X = {u} 
with Y being the variety of integro-differential algebras over K. 

Some identities following from S describe basic interactions between operations 
in the pure Baxter axiom (8), multiplicativity of the evaluation (5), the identities 

e 2 = e, «9e=0, eJ = 0, /(■/)* = (■/)/*, (24) 

and the variant (7) of the differential Baxter axiom connecting all three operations. 

We need to introduce some notational conventions. We use f,g for coefficients 
in and V for terms in i%(J^" U {«}). As for differential polynomials, we write 
u n for the nth derivative of u. Moreover, we write 

V(0)for E (V) and u(0) a for J^it l -(0) a ' J 

i=0 

where a is a multi-index. 

As a first step towards canonical forms, we describe below a system of terms 
that is sufficient for representing every integro-differential polynomial (albeit not 
uniquely as we shall see presently). 

Lemma 36. Every polynomial in JP{u} can be represented by a finite K-linear com- 
bination of terms of the form 

/«(0)Vj/i«W J- ••//„«*, (25) 

where f,f\ , . . . ,f n £ and each multi-index as well as n may be zero. 

Proof. The proof is done by induction on the structure of terms, using the above 
identities (8), (5), (20) and (24) of integro-differential algebras. □ 

With the aid of the previous lemma we can determine the constants of ,^{u}. 

Proposition 37. Every constant in ,^{u\ is represented as a finite sum £ a c a u{Q) a 
with constants c a in & ' . 

Proof. By the identity jd = 1 — e, a term V represents a constant in &{u} if and 
only if e(V) = V. Since V is congruent to a finite sum of terms of the form (25) and 
since Im(E) = c <o , the identities for e imply that V is congruent to a finite sum of 
terms of the form c a u(0) a . □ 

The above representation (25) of the integro-differential polynomials is not 
unique since for example when trying to integrate differential polynomials by using 
integration by parts, terms like 
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Ifvl and fu-Sfu-f(0)u(0) 
are equivalent. It becomes even more tedious to decide that, for instance, 

2x«(0)( 3 > 1 V 1 ' 3 ' ' 4 ) J(2* 3 + 3jc) u^J(x + 2) u^ 

and 

4jc«(0)( 3 ' 1 ) M ( 1 ' 3 ' ' 4 )J^ M ( 1 ' 2 - 3 )Jjc«( 2 )+6x«(0)( 3 > 1 ) M ( 1 > 3 A4)J xm (i,2,3)j^ + 2)m (2) 
+ 12x M (0)( 3 ' 1 ) M ( 1 - 3A4 )Jx«( 1 - 2 - 3 )Jm( 2 ) 
represent the same polynomial. In general, the following identity holds: 
Lemma 38. We have 

$Vu^u k+l = ^ (yu** 1 JV'«A +1 -V(0)u k {0) h+l ) (26) 

where k, fa > 0. 

Proof. Using (7) and the Leibniz rule, the left-hand side becomes 

S(v4*)(u k y vu£ +i - jv«* +1 - &JV«£*«* +] - v(o) u k {Q) h+ \ 

and the equation follows by collecting the JVk£* m^+i terms. □ 

The important point to note here is that if the highest derivative in the differential 
monomial ifi of order k+l appears linearly, then the term jfu^ is congruent to a 
sum of terms involving differential monomials of order at most k. This observation 
leads us to the following classification of monomials; confer also [10, 35]. 

Definition 39. A differential monomial u$ is called quasiconstant if j3 = 0, quasi- 
linear if j3 ^ and the highest derivative appears linearly; otherwise it is called 
functional. An integro-differential monomial (25) is classified according to its outer 
differential monomial uP, and its order is defined to be that of uP . 

Proposition 40. Every polynomial in JF{u} can be represented by a K-linear com- 
bination of terms of the form 

/«(0)Vj/i«»J-J/„«*, (27) 

where f,f\,...,f n G the multi-indices 0£,/3 as well as n may be zero and the 
lift , it 1 " are functional. 

Proof. By Lemma 36 we can represent every polynomial in ,^{u} as a /f-linear 
combination of terms of the form 

/«(0)Vj/mW J- ••//„«* (28) 
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where the multi-indices and « can also be zero. Let us first prove by induction on 
depth that every term can be written as in (28) but with nonzero multi-indices jk- 
The base case « = 1 is trivial since J f\ can be pulled to the front. For the induction 
step we proceed from right to left, using the identity 

j'/jv J/-JV JVJ7 

implied by the pure Baxter axiom (8). 

For proving that every multi-index Yk m (28) can be made functional, we use 
noetherian induction with respect to the preorder on J = jfiu^' J ■ ■ ■ SfnU?" that first 
compares depth and then the order of u^ x . One readily checks that the left-hand side 
of (26) is greater than the right-hand side with respect to this preorder, provided that 
V is of this form. 

Applying Lemma 38 inductively, a term J/im 1 ' 1 is transformed to a sum of terms 
involving only integral terms with functional differential monomials, and the base 
case n = 1 follows. As induction hypothesis, we assume that all terms that are 
smaller than J = fu(0) a uP Jftu* 1 J • • • J can be written as a sum of terms in- 
volving only functional monomials. Since J/2«^ J • • • J/nW^" is smaller than J, it 
can be written as sum of terms involving only functional monomials; we may thus 
assume that u^,... ,uY" are all functional. Since yi is nonzero, we are left with the 
case when is quasilinear. Applying again Lemma 38 inductively, we can replace 
u7 l in J by a sum of terms involving only integral terms with functional differential 
monomials. The induction step follows then by the linearity of J. □ 

For implementing the integro-differential polynomials in TH30REMV we use the 
functor hierarchy described in Sect. 2. The multi-index representation u$ for terms 
of the form (19) is realized by the monoid N* of natural tuples with finitely many 
nonzero entries, generated by a functor named TuplesMonoid. The nested integrals 
J/im 1 ' 1 J • • • jf n uY" are represented as lists of pairs of the form (fk,Yk)> with A € 
and % G N*. The terms of the form (25) are then constructed via a cartesian product 
of monoids as follows: 

Definition [ "Term Monoid for IDP", any[:F, N] , 
TermMonoid[:F, N] = TuplesMonoid[N] x TuplesMonoid[N] x TuplesMonoid [T x TuplesMonoid [N] ] ] 

Using this construction, the integro-differential polynomials are built up by the 
functor FreeModule[J^,B] that constructs the J^-module with basis B. It is instan- 
tiated with & being a given integro-differential algebra and B the term monoid just 
described. We will equip this domain with the operations defined as below, using a 
functor named IntDif f Pol[J^,K]. Later in this section we will present some sam- 
ple computations. 

Canonical Forms for Integro-Differential Polynomials. It is clear that ^-linear 
combinations of terms of the form (27) are still not canonical forms for the integro- 
differential polynomials since by the linearity of the integral, terms like 



f$(g+ h ) u and fh u +fP m 
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or terms like 

/ jXgu and Xf jgu 

with f,g,h £ & and X € K represent the same polynomial. To solve this problem, 
we can consider terms of the form (27) modulo these identities coming from lin- 
earity in the "coefficient" / and the integral, in analogy to the ideal [ introduced 
in Sect. 4 for J£" # [<9, J*]. Confer also [39], where the tensor product is employed for 
constructing free objects in differential Rota-Baxter algebras. In the following, we 
assume for simplicity that & is an ordinary integro-differential algebra. 

More precisely, let & denote the set of terms of the form (27) and consider the 
free Zf-vector space generated by 01. We identify terms 

fu(o) a uPjf lU *S---$f n » r " 

with the corresponding basis elements in this vector space. Then we factor out the 
subspace generated by the following identities (analogous to the construction of the 
tensor product): 

./'.'J./,»M ---.IX/A - .AikM -".!'./^ 

= fujft >r .\' r itw + .trif. «*;••■ 

fUjfm* J • • • JW*)"* J • • • Ir> r " = A//.' J./; J ••• • • • }/„„■■ 

Here 17 is short for u(0) a u$ , and there are actually two more identities of the same 
type for ensuring /^-linearity in /. We write \S%\ for this quotient space and denote 
the corresponding equivalence classes by 

Lf«(0)Vj/i«»J-J/„«*]. (29) 
By construction, the quotient module now respects the linearity relations 

L/i/j/^;... ja* +/*>*;• ••;/,«*] 

= /rj/, «« S-SM- S^] + [fU j flU n J . . . J/, J . . . ;/„«*] 
»' J • • • J(AA)«* J • • • J/«» r "] = ^ ./'''.I'/i JfkU % ! ■ ■ ■ //„«*]. 

together with the ones for linearity in /. 

As for the tensor product, we have canonical forms for the factor space by ex- 
panding the "coefficient" / and all the in (29) with respect to a K-basis 38 for 
assuming S3 contains 1. Then every polynomial can be written as a ^-linear 
combination of terms of the form 

M0) a M /3 jM ri /--->,,M rn , (30) 

where b,b\,...,b n € S§ with the condition on multi-indices as in Proposition 40. 

To show that terms of the form (30) are canonical forms for the integro-differential 
polynomials, we endow the quotient space with an integro-differential structure 
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and invoke Proposition 34. For this we define the operations on the generators (29) 
and check that they respect the above linearity relations on [M\ . 

First, we define a multiplication on [M]. Let &q C ffl denote the /if-subspace 
generated by integral terms \f\ u n J • ■ ■ $f n u7", including 1 S 29, as the case « = 0. 
Clearly, the nesting operation (23) can be defined in a completely analogous manner 
on such integral terms (the only difference being that we have now derivatives of 
the indeterminate). Since it is clearly ^-linear, it induces an operation U : [Mq] x 
[£%o] — > [&o]- The next step is to define the shuffle product on 2?q just as in (22), 
again with obvious modifications. Passing to the quotient yields the shuffle product 
• : [£%o] x [2?q] — > [2?q\. This product is finally extending to a multiplication on all 
of \3£\ by setting 

\fu(0) a uh][fu(0fuh~\ = \ffu{Q) a+& u^{J-J)] 
where J and / range over Let us compute an example: 

MultlDP^IO]* 1 '^ 2 ) "J" (3x) u* 1 ' 1 ' "_f"x 2 u<°' 2 >, 3u[0] <2 ' 3> u< 3 ' 1 > " J"" x u< 1 ' ' l >] 

-18 u[0] <3 ' 3> u* 5 ' 1 * Jxu"' 11 JV°' 2> Jx 2 u<°' 2 > - 27 u[0] <3 ' 3> u* 5 ' 11 J x u (1 ' x > J x 2 u <0 ' 2 > J u<°' 2 > + 
27 xu[0] <3 ' 3> u< 6 ' 2> Jxu' 1 ' 1 * Jx 2 u<°' 2 > - i xu[0] <3 ' 3> u< 7 ' x > J u<°' 2 > J x 2 u<°' 2 > - 

27u[0]< 4 ' 4 ' u' 5 ' 1 * Jxu* 1 ' 1 ' JVu<°' 2 > + ^u[0] <5 ' 3 > u* 5 ' 1 * JV°' 2 > JVu<°' 2 > 

Since the multiplication on & and the shuffle product are commutative, associative, 
and distributive over addition, the multiplication on [M] is well-defined and gives 
the structure of a commutative TT-algebra. 
The definition of a derivation d on this algebra is straightforward, using the fact 
that it should respect /f-linearity and the Leibniz rule (treating also the u(Q) a as 
constants), that it should restrict to the derivation on differential polynomials (which 
in turn restricts to the derivation on j£"), and finally that it should also satisfy the 
section axiom (3). Here is a sample computation: 

DiffIDP[u[0] <1 > u< 2 '!> »r.. xu <l.°.2> "J"" (3x 2 ) u<°' 2 > +3u[0] <2 ' 3> u< 3 ' 2 > "J"" (2 x 3 + 4 x) u< 2 ''>] 

6u[0] a> u' 1 ' 2 ' Jxu* 1 ' ' 2 ' JVu <0 ' 2 > + 3u[0] <1> u <2 '°' 1 > Jxu' 1 ' ' 21 Jx 2 u <0 ' 2> + 

3xu[0] (1> ui 3 ' 1 ' 2 * Jx 2 u<°' 2 > + 18u[0] <2 ' 3> u <2 ' 3 > JVu' 2 ' 1 * + (24 x + 6 x 3 ) u [ 0] <2 ' 3> u< 5 ' 3 > + 
12u[0] !2 ' 3> u< 3 ' 1 ' 1 > Jx 3 u< 2 ' 1> + 8ku[0] <2,S> u<«' 1 ' 1 > -3u[0]< 5 ' 3 ' u< 2 ' 3 > - 2 u [0] <5 ' 3> u* 3 ' 1 ' '> 

Using the ^-linearity of this derivation, one verifies immediately that it is well- 
defined. From the definition it is clear that ^-linear combinations of generators of 
the form [w(0) a ] are constants for d, and one can also check that all constants are 
actually of this form. 

Finally, we define a /f-linear integral on the differential /^-algebra ( ,d). Since 
we have to distinguish three different types of integrals, here and subsequently we 

will use the following notation: the usual big integral sign / for the integration to 

be defined, the small integral sign J for the elements of & as we have used it before, 
and /jf for the integral on 
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The definition of the integral on is recursive, first by depth and then by order 
of u$ , following the classification of monomials from Definition 39. In the base case 
of zero depth and order, we put 

/[/«(on=[/^/]Ko) a ]. (3D 

Turning to quasiconstant monomials, we use the following definition (which actu- 
ally includes the base case when / = 1): 

J [fu{o) a j] = nofiym Horj(yf)j'}. (32) 

In the quasilinear case we write the generators in form 

\fu{Q) a Vul k u k+l J] with V = Mq° • • • u^zl 

and construct the integral via (26). Writing s = jSj + 1, we have u^u^i = (uV)'/s, 
so we can define 

J [fu{Q) a V(u{)'j] = [fu{Q) a VulJ]~[u(0) a ] J [/V/]K]-/(0)[«(0r+^1, (33) 

where we write f(0) for e(/) and J is 1 for 7=1 and zero otherwise. In the func- 
tional case, we set 

j[fu(0) a uh] = [u(0) a Jfuh], (34) 

so here we can just let the integral sign slip into the equivalence class. One may 
check that the integral is well-defined in all the cases by an easy induction proof, 
using ^-linearity of the integral, the evaluation on J?, and the derivation on [&\ . 

Here is a small example of an integral computed in the quasiconstant case (note 
that IntlDP corresponds to the big integral and "J" to J in our notation): 

IntIDP[u[0]< 1> "J"xu< 1 '°' 2 > "J" (x 2 + 2) u<^ 2 >] 

2xu[0] u> Jxu< 1 '°' 2 > JV 1 ' 2 * +xu[0] a> Jxu <1 '°' 2 > Jx 2 u< 1 ' 2 > - 2 u [ ] <x> J x 2 u' 1 ' °' 2> J u< 1 ' 2 > - 
u[0] <x> JVu"' ' 2 ' JVu (1 ' 2 > 

The next example computes an integral in the quasilinear case: 

IntIDP[u[0] < 3 ' 2> u< 2 ^> "J"xu< 1 '°' 2 > "J" x 2 u^' 1 )] 

-u[0] <5 ' 2> Jxu<i' ' 2 > - -x 2 u[0] <3 ' 2> u< 2 > Jxu<"'°' 2 > + -x 2 u[0] <3 ' 2> u< 5 > Jxu' 1 ' ' 2 ) - 
6 6 6 

iu[0] <5 ' 2> u< 3 > [xu' 1 ' ' 2 ' 
6 

Note that all differential monomials within integrals are functional again, as it must 
be by our definition of [&] . 

By construction the integral defined above is a section of the derivation on . 
So for showing that \M] is an integro-differential algebra with operations, it remains 



46 M. Rosenkranz, G. Regensburger, L. Tec, and B. Buchberger 

only to prove the differential Baxter axiom (4). Equivalently, we can show that the 
evaluation 




is multiplicative by Corollary 4. 

Recall that the algebra of constants ^ in (\i%\,d) consists of /if-linear combina- 
tions of generators of the form [m(0)"]. By a short induction proof, we see that 

J[u{Q) a ][R] = [u{Q) a ]j[R]. (35) 

Hence the integral is homogeneous over the constants. 

For showing that the evaluation is multiplicative, we first reassure ourselves that 
it operates in the expected way on integro-differential monomials. 

Lemma 41. We have 

E\fu(0) a uh]=f(Q)[u(Q) a+ hl 
where J is 1 for J — 1 and zero otherwise as in (33). 

Proof. Note that e is ^-linear by (35), so we can omit the factor u(0) a . Assume 
first j3 = 0. Then by the quasiconstant case (32) of the definition of the integral, we 
have 

■ \fj] = \fj\ J if j}' = \fA [{yf)A + J i(yf)j'} - J i/n 

which by J>/' = / - /(0) gives 

f(0)[J]-f(0)j[J]' = f(0)[J\ 

because 

J [J]' = [J] for J^l 

by the functional case (34) and zero for J = 1 . If j3 ^ is of order k, we write 
ifi = Vul with s 7^ 0, and we compute 

e \fuh) = \fVu\J] - J [fVu s k J]' = /(0) [u(Q)Pf\ 
by the quasilinear case (33) and the Leibniz rule. □ 

Theorem 42. With the operations defined as above, ( \£%\ , d , J ) has the structure of 
an integro-differential algebra. 

Proof. As mentioned above, it suffices to prove that e is multiplicative, and we need 
only do this on the generators. Again omitting the u(Q) a , we have to check that 
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E [fuh] [fjf\ = E [ff U P + P (J ■J)) = E \fuh] • E \fjj\ . 

The case J = J = 1 follows directly from Lemma 41 and the multiplicativity of e 
in & . Otherwise the shuffle product J ■ J is a sum of integral terms, each of them 
unequal one. Using again Lemma 41 and the linearity of e, the evaluation of this 
sum vanishes, as does e[/m^7] • e[/m^7]. □ 

Since is an integro-differential algebra, we can conclude by Proposition 40 
and Proposition 34 that \£%\ leads to canonical forms for integro-differential poly- 
nomials, up to the linearity relations: After a choice of basis, terms of the form (30) 
constitute a system of canonical forms for JP{u}. In the TH30REMV implementa- 
tion, we actually compute in and do basis expansions only for deciding equality. 



7 From Rewriting to Parametrized Grobner Bases 

Equipped with the integro-differential polynomials, we can now tackle the task of 
proving the convergence of the reduction rules in Table 1 . As explained in Sect. 4, 
we will use the Diamond Lemma (Theorem 22) for this purpose. First of all we must 
therefore construct a noetherian monoid ordering > on JP(d, J) that is compatible 
with the reduction rules. In fact, there is a lot of freedom in defining such a >. It is 
sufficient to put d > / for all / € & and extend this to words by the graded lexi- 
cographic construction. The resulting partial ordering is clearly noetherian (since it 
is on the generators) and compatible with the monoid structure (by its grading). It 
is also compatible with the rewrite system because all rules reduce the word length 
except for the Leibniz rule, which is compatible because d > f. 

Thus it remains to prove that all ambiguities of Table 1 are resolvable, and we 
have to compute the corresponding S-polynomials and reduce them to zero. On the 
face of it, there are of course infinitely many of these, suitably parametrized by 
f,g G & and <p, \y € <P. For example, let us look at the minimal fork generated by 
JmJvJ. In this case, the rule J/ J may be applied either with / = u or with / = v 
yielding the reductions 

/ \ 

(J.«)JvJ-J(J.«)vJ J"(J-v)J-J«J(J-v) 

with the S-polynomial p = (J ■ u) JvJ - J (J • u)vj - J« (J • v) J + JmJ (J • v). But 
actually we should not call p an S-polynomial since it represents infinitely many: 
one for each choice of u, v € 

How should one handle this infinitude of S-polynomials? The problem is that 
for reducing S-polynomials like p one needs not only the relations embodied in the 
reduction of Table 1 but also properties of the operations d , J : — >• & acting on 
u,v G & . Since these computations can soon become unwieldy, one should prefer a 
method that can be automated. There are two options that may be pursued: 
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• Either one retreats to the viewpoint of rewriting, thinking of Table 1 as a two- 
level rewrite system. On the upper level, it applies the nine parametrized rules 
with f.g £ J? and (p,yr E <P being arbitrary expressions. After each such step, 
however, there are additional reductions on the lower level for applying the prop- 
erties of <9, J: & — > & on these expressions. Using a custom-tailored reduction 
system for the lower level, this approach was used in the old implementation for 
generated an automated confluence proof [70]. 

• Or one views an S-polynomial like p nevertheless as a single element, not in 
J^"(<9, J) but in ^(d, J) with JF = JF{m,v}. With this approach, one remains 
within the paradigm of parametrized Grobner bases, and the interlocked opera- 
tion of the two levels of reduction is clarified in a very coherent way: The upper 
level is driven by the canonical simplifier on J^[<9, J], the lower level by that on 
^{u,v}. 

It is the second approach that we will explain in what follows. 

Using JF(d, J") instead of ^(d, J") takes care of the parameters f,g £ & but 
then there are also the characters (p, yf £ <P. The proper solution to this problem 
would be to use a refined version of integro-differential polynomials that starts from 
a whole family (J 9 ,)<pe* of integrals instead of the single integral J, thus leading 
to a corresponding family of evaluations u((p) instead of the single evaluation m(0). 
We plan to pursue this approach in a forthcoming paper. For our present purposes, 
we can take either of the following positions: 

• The characters <p, yf may range over an infinite set <P, but they are harmless since 
unlike the /,g£# they do not come with any operations (whose properties must 
be accounted for by an additional level of reduction). In this case, Table 1 is still 
an infinitely generated ideal in £F{d, J), and we have to reduce infinitely many 
S-polynomials. But the ambiguities involving characters are all of a very simple 
nature, and their reduction of their S-polynomials is straightforward. 

• Alternatively, we may restrict ourselves to a finite set of characters (as in most 
applications!) so that Table 1 actually describes a finitely generated ideal in 
J?(<9, J), and we need only consider finitely many S-polynomials. 

The second alternative is somewhat inelegant due to the proliferation of instances 
for rules like (pyr — > yf. In our implementation, we have thus followed the first alter- 
native with a straightforward treatment of parametrization in <p , yf but we will ignore 
this issue in what follows. 

We can now use the new TH30REMV implementation for checking that the nine 
rules in Table 1 form a Grobner basis in J£"(<9, J). As explained before, we use 
the Diamond Lemma for this purpose (note that the noetherian monoid ordering > 
applies also to &(d, J) except that we have now just two generators u,v £ = 
J£"{m,v} instead of all / £ JF). Hence it remains to check that all S-polynomials 
reduce to zero. We realize this by using the appropriate functor hierarchy, as follows. 
We first build up the algebra of the integro-differential polynomials having, in turn, 
integro-differential polynomials as coefficients, via the functor 

IntDiffPolys[IntDiffPolys[^,K],K] 
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and we denote the resulting domain by P. Then we consider an instance of the func- 
tor constructing the integro-differential operators over P. Finally, the computations 
are carried out over the algebra created by the GroebnerExtension functor taking 
the latter instance as input domain, that allows to perform polynomial reduction, 
S-polynomials and the Grtibner basis procedure. 

Of course, the S-polynomials are generated automatically, but as a concrete ex- 
ample we check the minimal fork considered above: 

ReducePol[ ( ("_f" . u (1 >) "J"" . v« "J"" - "J"" ("J" . u (1 >) v« "J"") - 

( „r„ u <i> („;„.„<!>) „;„ _.;.„<» „;„ (-/..,«))] 

o 

As it turns out, there are 17 nontrivial S-polynomials, and they all reduce to zero. 
This leads us finally to the desired convergence result for JF[d, J"]. 

Theorem 43. The system of Table 1 represents a noncommutative Grobner basis in 
&(d, J) for any graded lexicographic ordering satisfying d > f for all f E 

Proof By the Diamond Lemma we must show that all S-polynomials p E J^"(<9, J) 
reduce to zero. Since they may contain at most two parameters f,g E let us write 
them as p(f,g). But we have just observed that the corresponding S-polynomials 
p(u,v) E J?(<9, J) with & = J?{u,v} reduce to zero. Using the substitution homo- 
morphism 

q>: (u,v)^(f,g), 

lifted to ^[d,j] — > JP[d,j] in the obvious way, we see that p(f,g) = (pp(u,v) 
reduces to zero as well. □ 

From the conclusion of the Diamond Lemma, we can now infer that Table 1 indeed 
establishes a canonical simplifier for JF[d,j]. 



8 Conclusion 

The algebraic treatment of boundary problems is a new development in Symbolic 
Analysis that takes its starting point in differential algebra and enriches its structures 
by introducing an explicit notion of integration and boundary evaluations. Recall the 
three basic tools that we have introduced for this purpose: 

• The category of integro-differential algebras (^,d,j) for providing a suitable 
notion of "functions". (As explained in Sect. 2, here we do not think of categories 
and functors in the sense of Eilenberg and Maclane — this is also possible and 
highly interesting but must be deferred to another paper.) 

• The functor creating the associated integro-differential operators &[d, j] as a 
convenient language for expressing boundary problems (differential operators, 
boundary operators) and their solving Green's operators (integral operators). 
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• The functor creating the associated integro-differential polynomials ^{u}, which 
describe the extension of an integro-differential algebra by a generic function u. 

In each of these three cases, the differential algebra counterpart (i.e. without the 
"integro-") is well-known, and it appears as a rather simple substructure in the full 
structure. For example, the differential polynomials ,^{u] over a differential alge- 
bra (J?, d) are simple to construct since the Leibniz rule effectively flattens out com- 
pound terms. This is in stark contrast to an integro-differential algebra (J^,<?, J), 
where the Baxter rule forces the presence of nested integrals for ensuring closure 
under integration. 

The interplay between these three basic tools is illustrated in a new confluence 
proof: For an arbitrary integro-differential algebra (J^ - , d, J), the rewrite system for 
the integro-differential operators & [d , j] is shown to be a noncommutative Grobner 
basis by the aid of the integro-differential polynomials ^{u,v}. Having a confluent 
rewrite system leads to a canonical simplifier, which is crucial for the algorithmic 
treatment as expounded in Sect. 2. 

Regarding our overall mission — the algebraic treatment of boundary problems 
and integral operators — we have only scratched the surface, and much is left to be 
done. We have given a brief overview of solving, multiplying and factoring bound- 
ary problems in Sect. 5. But the real challenge lies ahead, namely how to extend our 
framework to: 

• Linear Boundary Problems for LPDEs: As mentioned at the start of Sect. 5, the 
algebraic framework for multiplying and factoring boundary problems is set up 
to allow for LPDEs; see [66] for more details. But the problematic issue is how 
to design a suitable analog of &[d, J] for describing integral and boundary op- 
erators (again the differential operators are obvious). This involves more than 
replacing d by d/d x ,d/d y and J by J Q , j } because even the simplest Green's 
operators employ one additional feature: the transformation of variables, along 
with the corresponding interaction rules for differentiation (chain rule) and inte- 
gration (substitution rule); see [74] for some first steps in this direction. 

• Nonlinear Boundary Problems: A radically new approach is needed for that, so it 
seems appropriate to concentrate first on boundary problems for nonlinear ODEs 
and systems thereof. A natural starting point for such an investigation is the dif- 
ferential algebra setting, i.e. the theory of differential elimination [44, 1 1, 12]. By 
incorporating initial or boundary conditions, we can use explicit integral opera- 
tors on equations, in addition to the usual differential operators (prolongations). 
As a consequence, the natural objects of study would no longer be differential 
but integro-differential polynomials. 

We are well aware that such an approach will meet with many difficulties that 
will become manifest only as we progress. Nevertheless, we are confident that an 
algebraic — and indeed symbolic — treatment along these lines is possible. 
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